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Abstract. The confinement transition from the quark and gluon degrees 
of freedom appropriate in perturbation theory to the hadrons observed by 
real world experiments is poorly understood. In this strongly interacting 
transition regime we presently rely on models, which to varying degrees 
reflect possible scenarios for the QCD dynamics. Because of the absence 
of beam and target remnants, and the clean experimental conditions and 
high event rates, e + e _ annihilation to hadrons at the Z° provides a unique 
laboratory, both experimentally and theoretically, for the study of parton 
hadronization. This review discusses current theoretical understanding of 
the hadronization of partons, with particular emphasis on models of the non- 
perturbative phase, as implemented in Monte Carlo simulation programs. 
Experimental results at LEP and SLC are summarised and considered in the 
light of the models. Suggestions are given for further measurements which 
could help to produce more progress in understanding hadronization. 



Topical Review to be published in 

Journal of Physics G: Nuclear and Particle Physics 

Short title: Hadronization in Z° decay 

February 1, 2008 

|| On leave from Department of Physics and Astronomy, The University of Manchester, 
Manchester, M13 9PL, UK. 



2 



1. Introduction 

Hadronic systems produced in e + e~ annihilation have their origin in a uniquely 
simple quark- ant iquark state. While the standard model of particle physics 
provides a well tested description of the reaction e + e~— > Z°/7* — > qq, the 
subsequent production of observable hadrons is less well understood. A 
parton shower, described by perturbative quantum chromodynamics (pQCD), 
is normally invoked to describe the initial fragmentation phase. In the 
subsequent non-perturbative hadronization process, the partons become the 
hadrons which are experimentally observed. Multihadronic e + e~ annihilation 
events, with no beam or target fragments to confuse their experimental or 
theoretical interpretation, provide the most powerful system available for the 
study of the transition from the partons of perturbative QCD to the hadrons 
of the laboratory. 

Between 1989 and 1995, the Large Electron Positron collider (LEP) at 
CERN delivered an integrated luminosity of about 170 pb -1 per detector at and 
around the Z° peak, providing each of its four dedicated experiments, ALEPH, 
DELPHI, L3 and OPAL, with some 6 million e + e~ annihilation events, 70% 
of which were multihadronic. These events have enabled the experiments to 
conduct detailed studies of many aspects of parton hadronization. Over the 
same period, the SLAC Linear Collider (SLC) delivered a more modest number 
of events, initially to the MARK II detector and then to the SLD detector, 
which has now accumulated about 200k Z° events. Although its luminosity is 
lower than that of LEP, the SLC is able to provide highly polarized electron 
beams. While LEP has now entered phase 2, running at higher energies, 
analysis of the Z° peak data will continue for some time. The SLC, with 
polarized beams, is scheduled to run at the Z° for several more years. 

In this review, we look at the present understanding of the non-perturbative 
hadronization process in the context of the recent experimental results 
from LEP and SLC. A brief introduction to the electroweak theory of 
e + e~ annihilation is followed by a heuristic picture of the transition from 
perturbative partons to final-state hadrons. After this, the role of perturbative 
QCD in determining the structures and particle content of events is discussed. 
Models of the non-perturbative hadronization phase are then covered - 
independent fragmentation, string fragmentation and cluster fragmentation 
are discussed in detail and a comparison is given of the main features of the 
different models. The LEP and SLC machines and their associated detectors 
are then described, with emphasis on the features relevant to measurements 
of individual hadrons from Z° decay. Next the experimental data on single 
inclusive identified particles are summarised, and the results interpreted in 
the light of the different models. Results on spin phenomena, such as baryon 



3 



polarization and vector meson spin alignment are then covered, and this is 
followed by discussion of results on correlations, including Bose Einstein effects, 
strange particle rapidity and angular correlations, and intermittency; where 
possible the results are interpreted within the context of theory and models. 
Differences between quark and gluon initiated jets are then dealt with, and the 
review ends with a look forward to what may still be learned from the available 
LEP and SLC data, and the forthcoming SLC data. 



1.1. Electroweak aspects of e + e annihilation 

The (initial) numbers, directions, polarizations and flavours of the quarks 
produced in e + e~ annihilation are determined by the electroweak couplings of 
the exchanged vector bosons. The full, tree-level expression for the differential 
cross section for e + e~^ qq in the centre of momentum (CoM) frame is given 
as a function of quark's polar angle, 6*, measured relative to the electron beam 
direction, and total CoM energy squared, s, by: 

^ (S ' r) = ^£^ q £ XiWxMtivi + + x 




(2 - 7>W)Bj'(l - Pi) sin 2 e*+(v™ + V®A%) (1 + cos 2 0*) 



+ (V^Al j + Vi 2) ) (V^A'l + V® )2 cos 6* j (1) 
where, 

<p(l)_l nn - Aij 



v\v j { +$a\4 w (s-M 2 ) + iM,r 



The summation is over all exchanged bosons, % = 7*,Z°, . . ., with mass and 
width Mi and Tj, and vector and axial couplings v\ and a\. For example 
v] = Q t , a] = 0, vf = (T f 3L - 2Q iS 2 w )/2s w c w and af = Tf L /2s w c w 
etc, with Q{ the fermion's electric charge, normalized to that of the electron, 
Tf L the fermion's third component of weak, SU(2) L , isospin and sw(c\y) the 
sine(cosine) of the Weinberg angle. The number of colours is N c = 3 and a cm 
is the electromagnetic fine structure constant. The dependence on the initial 
lepton and final quark longitudinal polarizations are via pf, where pf = +1 
for spin along the CoM direction of travel (helicity basis). In the absence of 
transverse beam polarization there is no azimuthal angular dependence. 

At v^i = Mi the Z° exchange term dominates and (1) simplifies 
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considerably. In this limit the partial width for Z° 

N c a em M z 



qq is given by: 



qq 



6 



^ q (3 - P 2 q )< + 2fl& 



(3) 



This implies relative hadronic branching ratios of ~ 17% to each up-type quark 
and ~ 22% to each down-type quark at the Z°; the relative fraction of b quarks 
peaks on resonance. Figure p] shows the relative fractions of light (u,d,s), charm 
and bottom quarks as a function of \fs obtained using the full expression (1). 
The total hadronic branching fraction is ~ 70%. The beam polarization has 
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Figure 1. The relative branching fractions of light (u,d,s), charm and 
bottom quarks in e + e~ annihilation as a function of y/s for an unpolarized 
(solid lines) and 100% right polarized (dashed lines) electron beam. 



little influence on the relative rates of produced quarks, although the total rate 
is proportional to + a^)V^ — 2v e a e V^ . However only a single polarized 
beam can significantly alter the quark polar angle distribution, as can be seen 
from the forward-backward asymmetry (for pure Z° exchange and /3 q = 1): 

(-cos**) = V^A C -VP 2COS0* 
g(cos0*) + g(-cos0*) <pP-vPA e q l + cos 2 ^ {) 

Since A e is only w 0.16 (Ad ~ 0.94 and A u w 0.69) then, given polarized 
beams, one has a powerful statistical way to identify separately quark and 
antiquark jets based on their direction. The quarks are naturally produced 
highly polarized: 

mi* = + 1 ) ~ = - 1 ) - A i 4,(l-^)2cosg* ^ 

m (Pf = + 1 ) + f(Pf = -!) " l + cos 2 ^ + A c A q 2cos^ ~ q y> 



5 



This result is for unpolarized lepton beams, pure Z° exchange and /3 q = 1; the 
full expression depends little on the lepton beam polarization. 

For QCD studies at the Z° [|I| the effect of initial-state electromagnetic 
radiation (ISR) is of rather minor significance, unlike in the case of electroweak 
studies P, 0| . The primary effect of ISR is to lower the effective y/s, leading 
to a distortion of the Breit-Wigner lineshape given by (1). However the basic 
quantities of relevance from a QCD perspective, quark flavour mix, polarization 
etc, are only rather weak functions of y/s. Further, because of the Z° resonance 
the cross section falls rapidly when ISR occurs, mitigating against its effects 
(even so the peak cross section falls by ~ 25% compared to (1)), in contrast 
to the situation at LEP 2 and higher energies where ISR is very important @] . 
The effects of final-state radiation, particularly of gluons but also of photons, 
form one theme of this review. The effect of final-state radiation on the total 
hadronic cross section can be summarized in the multiplicative factor P, ||: 



where a s is the strong coupling constant (g^/Air) and Cp = 4/3 is a measure 
of the quark-gluon coupling strength. Flavour (mass) dependent effects and 
electroweak radiative corrections are small — both occur at the 1% level, and 



1.2. A picture of the parton to hadron transition 

The description of a multihadronic event commences by specifying a set of 
primary partons (q/q/g) distributed according to an exact Z° decay matrix 
element, such as (1). It is customary to identify three basic stages in the 
transition of these partons into detected hadrons ||. First a parton (or, 
equivalently, dipole) cascade, formulated according to pQCD, evolves the 
primary partons from the hard scattering scale Q ~ Mi into secondary partons 
at a (fixed) cut-off scale Qo ~ 1 GeV. It is during the calculable stages, hard 
subprocess and shower, that the event's global features are determined: energy 
dependences, event shapes, multiplicity etc. In a second stage, carried out at 
the fixed, low virtuality scale Q , a model is employed to convert the secondary 
partons into hadrons. The second stage is essentially energy (Q) independent, 
up to power corrections, and assumed to be local in nature. Finally unstable 
primary particles are decayed into stable hadrons and leptons according to 
decay tables [[jj. Schematically the fragmentation function is given by: 




(6) 



need not concern us here [[|. 



D£(x, Q 2 ) = (pQCD evolution: Q 2 -> Qjj) <g> (model:a -> H) 



(7) 



<g> (tables:H -> h,h', ..) 
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Here, D£(x, Q 2 ) is the probability to find a hadron of type h carrying a fraction 
x of the parton's momentum, in a jet initiated by the parton a, whose maximum 
virtuality is Q. The H represent possible intermediate hadrons. 




Way ^2-200 fin 7 * / ' Z° t sep « (fc^i?)" 1 • t had 



Figure 2. A schematic diagram of the spatial evolution of a hadronic Z° 
decay indicating the relevant time scales and distances associated with gluon 
bremsstrahlung and hadron formation. 

We consider this process in terms of its space-time structure: the discussion 
is based on |J and illustrated in figure 0. The hard subprocess, e + e~— * qq, 
may be viewed as the production of a highly virtual photon, or a real Z°, 
which impulsively kicks a qq pair out of the vacuum; the time scale is short, 
^annih ~ 1 / ~ 10~ 2 fm. In this non-adiabatic process the quarks shake 
off most of their cloud of virtual particles, so that any structure they have is 
on a scale below \j\fQ 1 . That is they behave as bare (more properly half- 
dressed) colour charges until the gluon field has had time to regenerate out 
to a typical hadron size R m 1 fm. Allowing for the boost, this takes a 
time thad ~ Q/th x R ph QR 2 ~ 10 2 fm, where the second approximation 
is appropriate to light hadrons. 

The fact that thad t annih raises the issue of how charges are conserved 
over the space-like separated distances involved. Of course the accelerated 
quarks will radiate gluons and here two new time scales are relevant. First, 
from the off-shellness of the quark prior to emission, the formation time of 
a real gluon of 3-momentum (=energy) k is tf orm ~ kjk\ (the uncertainty 
relation gives the proper lifetime of the virtual state as the reciprocal of its off- 
shellness, during which time it will travel a 4-distance q^/q 2 in the laboratory 
frame, where g M is the virtual state's 4-momentum). Second, for the gluon 
to reach a transverse separation of R and become independent of the quark 
takes a time t sep = (k±R) ■ tf orm , whilst the hadronization time may be written 
£had(~ kR 2 ) = (k ± R) 2 ■ £ form . 

For this quark- gluon picture to make sense we require k± > R^ 1 so that 
iform < ^so P < ^had! this is natural as it implies that a s (/c 2 ) < 1, making 
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perturbation theory applicable. If k± < R^ 1 then we can say nothing. On 
the borderline are quanta with kj_ = R~ x ; these feel the strong interaction and 
are responsible for holding hadrons together. We distinguish these from the 
essentially free perturbative gluons by the name gluers. The first gluers form 
after only 1 fm, having k ~ k± « Br 1 . Thus gluers immediately form in the 
wake of the primary partons blanching their colour field and leaving two fast 
separating charge neutral systems. Thus the slowest hadrons form first, close 
to the interaction point, in what may be called an 'inside-out' pattern O. 

On quite general grounds the distribution of gluers in QCD can be 
estimated using: 



dA/l-lnnrs ~ Cf 



f ^s(kl) v ^dk 2 A 

Jk ± ~R IT 



k i 



dk dk 



Here the logarithmic k\ dependence reflects the fact that the coupling constant 
is dimensionless, the logarithmic k dependence follows from the gluon being 
massless and only the kernel function V, which is of order unity, depends 
on details of the quark-gluon vertex. The distribution is thus governed by 
longitudinal phase space and leads to hadron production on rapidity plateaux 
along the directions of the initial partons. 

Not all quanta are emitted at low k± and in a significant range, R- 1 < k± < 
k <C VQ 2 , a shower of perturbative gluons is possible. A hard gluon emitted 
at an angle 9 becomes a separate colour source at t sep = R/9(= (g±R) ■ t iorm ), 
again raising the issue of charge conservation. Fortunately a gluer emitted at 
the same angle 9 would appear just in the right place and at the right time 
to blanch the tail of the separating gluon's colour field. This gluon then acts 
as a new source of gluers, restricted such that k^9~ l R~ 1 . That is equivalent 
to a boosted (7 = 1/9) quark jet of a reduced scale Q = k± and with the 
substitution Cp Ca for the parton's charge in (||) ]l(J. It should be noted 
that gluers emitted at larger angles d > 9 have R/d < R/9 = t sep and should 
therefore be associated with the parent q+g system. For example at time R 
the first gluers emitted actually see ensembles of partons equivalent to the 
original qq pair. The coherence of soft emissions simply reflects colour charge 
conservation and is fundamental to a gauge theory such as QCD. 

Finally, the time scales for the decay of the primary particles are set by 
the reciprocal of their widths: for the strong resonances typical values of 
T ~ 1 — 100 MeV give tdecay ~ 200—2 fm. Clearly these scales are commensurate 
with those for primary hadron production so that any distinction between 
hadronization and resonance decays may be only semantics. This is illustrated 
in figure 13 which shows a string-inspired space-time picture (see section |3.2|) of 



two equivalent ways to form a 7r + 7r pair differing only in the time sequence 
of the string breaks. In the first figure the pions are directly produced whilst 
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in the second they appear to come from a p°; however observationally they are 
indistinguishable. The time scale for the weak decays of s, c and b hadrons is 
of order idecay ~ (M-w/mq) 2 k/mq (where k/rriQ is the hadron's 7 factor); this 
is ~ 10 12 fm for c and b and much longer for s, and so these decays may safely 
be treated as separate subprocesses. 



x 



Figure 3. Two possible views from the string model perspective of it + it~ 
production via an intermediate p°, illustrating the practical difficulty in 
separating strong resonance decays from the direct string fragmentation. 



1.3. Monte Carlo event generators 

A typical hadronic final state in Z° decay contains about 20 stable 
charged particles and 10 neutral hadrons, mainly 7r°'s. This is a complex, 
non-perturbative system beyond direct, first principles calculation at the 
present time. However Monte Carlo event generator programs || |TT| have 
been developed which provide remarkably accurate, detailed descriptions of 
complete hadronic events. These are based on a combination of judicious 
approximations to pQCD, discussed in section 0, and per/in-spirational models 
for hadronization, to be discussed more fully in section [| Here we list the main 
features of the major event generators used in Z° studies and refer the reader 
to later sections and to the literature for more complete descriptions. 



ARIADNE |12| implements a pQCD shower based on the dipole cascade 



model | 13fl , which is equivalent to a coherent parton shower. The 
evolving chain of dipoles corresponds naturally to a string and indeed 
the hadronization is performed using the JETSET implementation of the 



Lund string model |L4f . 



CO JETS (TJj is based on a virtuality-ordered parton shower, that takes 
no account of colour coherence. The resulting jets are hadronized 
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using a refined version of the Field-Feynman independent fragmenation 



model 16 



HERWIG [O] is based on a highly developed parton shower algorithm [JTR 
Bl, that automatically takes into account colour coherence, and a 



relatively simple cluster hadronization scheme |20| 



JETSET 



21] uses a virtuality-ordered parton shower with an imposed 

A 



angular ordering constraint to take into account colour coherence 
relatively sophisticated hadronization scheme based on the Lund string 
model is provided ||14|| . An option also exists which does not employ a 
parton shower but tries to use instead higher-order hard-subprocess matrix 
elements plus string hadronization. 



2. The role of perturbation theory in e + e annihilation to hadrons 

As stated in section |1.2| the bulk properties of hadronic events in Z° decay 
are established early in the fragmentation when virtualities are large and QCD 
perturbation theory is valid. It is an important issue to establish to what extent 
pQCD dominates and what contributions are made by non-perturbative effects. 



2.1. Matrix elements and showers 

Two basic approaches are used to calculate hadronic event properties. The 
use of fixed-order perturbation theory is justified by the smallness of the 
running strong coupling constant [^] at Mz, a s {Q 2 ) = 4n/ (3 \n(Q 2 / A 2 ) ps 
0.12. Known results include the total hadronic cross section to three loops, 
order-af ||. Complete one-loop, order-a;? calculations are available |23| for 



planar (i.e. 3-jet dominated) event shapes, including contributions due to 
the orientation of the event plane with respect to the beam direction |^4 | 



and to quark mass effects |[25|| . Partial order-a^ results are known for 4-jet 
distributions fStjfl . Tree-level calculations, up to order-af, of 5-jet distributions 



are also available [ 27 ] . However a complication arises in this approach because 
tree-level diagrams diverge whenever external partons become soft or collinear 
and related negative divergences arise in virtual (loop) diagrams (in addition 
to ultraviolet divergences). Fortunately, in sufficiently inclusive measurements 
such as the total hadronic cross section it is guaranteed that the two sets 



of divergences cancel ||28|| . Unfortunately in more exclusive quantities such as 
event shapes the cancellation is no longer complete and large logarithmic terms 
remain, generically of the form L = \ii(Q 2 /Qq). Since a s L is of order unity 
this can spoil the convergence of naive perturbation theory. 

In the second approach the perturbation series is rearranged in terms of 
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powers of a s L: S n a n (a s L) n + a s (Q) Y^n b n {&sL) n + • • -. The first infinite 
set of terms represent the leading logarithm approximation (LLA), then the 
(^-suppressed next-to-LLA and so on. Many of these terms have been 



identified |2I5j and summed using renormalization group techniques. They 
are conveniently expressed via Q 2 -dependent fragmentation functions whose 



evolution is controlled by Altarelli-Parisi type equations |fJU[ : 
2 dDZ(x,Q 2 ) _ /-i dz a s v b fx 2 \ 

Here, the so-called splitting functions, P£ h ,(z), may be thought of as giving 
the probability of finding parton b (and b') inside a and carrying a fraction z 
of its momentum. Monte Carlo event generators implement solutions to these 
equations as parton showers || [IT], [HJ]. More recently the resummed results 



for a number of event shapes have been calculated and combined with 
the fixed order approach |33fl . The Durham jet clustering algorithm [34[] for 
example was proposed to allow such a resummation for jet rates. 

2.2. Global event properties 

At very low little structure is present in hadronic events, beyond that due 
to the presence of hadronic resonances, and these can be described essentially 
by isotropic phase space This is because the momentum scales involved 



in the hadronization, typically 200 MeV, are comparable with those of the 
proto-jets. At higher energy, due to the preferred collinear nature of gluon 
radiation, the hadrons clearly form collimated pairs of back-to-back jets with 
angular distributions compatible with those of the primary qq pairs 



Around y/s ~ 30 GeV a significant fraction of events deviate from a simple 
linear configuration: this was interpreted as first evidence for gluons J37|. At 



the Z°, hard, non-collinear, gluon radiation is manifested in multi-jet event 
structures |38] . 



Event shape variables have been developed to be sensitive to the amount 
of acollinearity and hence the presence of hard gluon radiation 0. Since 
theoretical predictions are at parton level |p9| , hadronization corrections, 



derived from Monte Carlo models, must be applied before making comparisons 
to data. Away from phase space boundaries the corrections are modest, but not 
negligible (of order 10%), and order-a 2 calculations generally describe the data 
well | HI, particularly so when resummed [jfl ]. In the case of planar variables 



such as thrust, transverse momentum in the event plane, wide jet broadening 
etc, this allows a s to be measured, and good consistency is seen with other 
measurements [J42[| . Only a leading-order description of aplanar, four-parton, 
event shapes such as the D-parameter, transverse momentum out of the event 
plane p° ut , the four-jet rate etc, are presently available, although partial results 
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are available []26] in higher order. Here agreement between data and Monte 
Carlo is less satisfactory, with 30% discrepancies occurring for p^ ut and the 
4-jet rate [11], [13], |44|. Hopefully the situation will be improved by matching 
the showers to exact fixed-order results. This will prove important because 
of a reliance on Monte Carlo models in establishing variables to discriminate 
at higher energies between continuum QCD events and other physics signals, 
such as W + W~ pairs. This is in addition to their importance in QCD and in 
increasingly demanding electroweak measurements such as those involving jet 
charge |45] . 



Establishing agreement between the predictions of event generator 
programs and data for event shape variables has proved very important for 
confirming our ability to model QCD as a whole and has led to a first generation 
tuning of the Monte Carlo models [17]. A second generation of tunings, 
which also take into account data on identified particle production, are now 
becoming available [11, [13], 43]. This tuning is not a trivial exercise as it is 
often the case that a distribution depends in a complex way on a model's free 
parameters. However some specific sensitivities have been identified. The 3-jet 
rate is very sensitive to the QCD scale parameter A, indicating the dominance 
of the shower and pQCD. Event shapes and inclusive momentum spectra are 
sensitive to the shower cut-off Qq, and also to parameters controlling the 
generation of transverse momentum in the hadronization. In programs that 
lack a shower stage it is necessary to retune the hadronization model at each 
a/s since this must describe the whole fragmentation process including the 
perturbative y/s dependence found in a shower. 



2.3. Power corrections 

As noted above, hadronization corrections, which at present are inherently 
model dependent, need to be applied to partonic predictions before comparison 
can be made to hadronic event properties. Empirically the differences between 
quantities measured at the hadron level and the partonic predictions are found 
to be power behaved: 5X = Xh a d — A par ~ l/Q n - For example, in the case of 
massless quarks, 5-Rhad ~ l/<3 4 HH, wnere Rha,d is the usual ratio of hadronic 
to fi pair cross sections; in fact this result follows from the operator product 
expansion (OPE) j|9| and can be related to the value of the gluon condensate: 
<5-Rhad ~ (a s G.G) /Q A . In the case of average event shape variables, a 1/Q 
behaviour is found (the OPE does not apply in this case, due to the presence 
of multiple scales, and indeed no corresponding dimension-1 local operator 
exists). This behaviour may be established by going to the power-enhanced, 
low Q, region where a 1/Q n variation is easily distinguished from the slow 
logarithmic variation of the perturbative expressions. Such a 1/Q behaviour is 
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significant because a 2 s {Q) ~ 1 / 'Q\q =m ■ thus to be able to take advantage of a 
future 0(«f) prediction will require a better understanding of hadronization. 
Three sources of power corrections are known: instantons, infrared (IR) 



renormalons and ultraviolet renormalons P8| . Only IR renormalons are 



believed to be relevant to the issue of hadronization corrections; instanton 
effects are too highly power-suppressed. The coefficients of pQCD series in a s 
generically suffer from a factorial growth, ~ fi^nl, making the series formally 
divergent. A standard technique for dealing with this type of behaviour is 
to employ a Borel transformation f48|. A pole in the transform at n.27r//3o 
would correspond to an exp(— n.2ir/ f3 a s ) = (A/Q) n term in the Borel summed 
series: these poles are called renormalons. Unfortunately the 'residue' of the 
pole, equivalent to the coefficient of the (A/Q) n term, is not calculable by this 
purely mathematical technique. 

Recently a relationship has been suggested between the positions of the 
renormalons and the power behaviour of perturbative series in the presence 
of an IR 'regulator' The basic idea is that the (resummed) perturbative 
calculation probes regions of phase space involving low-virtuality partons where 
non-perturbative confinement effects should also be important. If these regions 
are isolated by introducing a 'cut-off' into the perturbative calculation then 
requiring that any cut-off dependence is compensated by the non-perturbative 
hadronization correction allows the power behaviour to be determined. Using 
this idea, together with a gluon mass ^TJ, a 1/Q correction was derived for 
the thrust T, C-parameter and longitudinal cross section ctl with all the 
coefficients proportional to a common scale: 5(T) oc — ACp/ir, 5(C) oc 6Cp 
and 5(<jl) Cp. Refinements of this calculation |52| have attempted to relate 
the common scale to a fixed a s , representing an effective measure of long-range 
confinement forces at an inclusive level. A dispersive approach indicates that 
the numerator of the 1/Q term may involve several coefficients, for example 
proportional to ln(Q/A). 

The subject of power corrections is an active theoretical area and 
still subject to dispute. However there exists the exciting prospect of a 
phenomenology of power corrections. This is particularly so if the corrections 
to various processes can be related and so shown to be universal |53| or could 
be developed for full distributions rather than average values. For example, 
according to the calculations above, (T + (2/37r)C) would have no leading 1/Q 
power correction, enabling a more accurate high-Q prediction |[I . 



It is interesting to compare these speculations with the results of a simple 
tube model |^2|, |53| calculation. The idea is that after hadronization a parton 
jet is equivalent to a tube of hadrons distributed uniformly in rapidity, y, along 
the jet axis with transverse mass /i/unit y and length Y. Thus the total energy 
and momentum of the jet are: E = fisinhY and P = fi(coshY — 1). In a two- 
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jet event the thrust is simply given by T = P/E, so that 5(T) = —2(a/Q; 
likewise S{C) = 3irfi/Q and 5{a L /a tot ) = (ir/2)fi/Q. Fitting to data, all these 
estimates are consistent with \i = 0.5 GeV. 



2. 4- Direct photons 

Direct photons, which are unaffected by hadronization, offer an alternative 
way to probe the early perturbative event structure ||55| . In these events, 
cross sections are not simply given in terms of the quark charges but reflect 
the competition in showers between q— >qg and q— >qy branchings in a way 
that is sensitive to the choice of evolution variable |56||. Here agreement with 
experiment has been less satisfactory |j57|, particularly for rather soft or very 



isolated photons; for example n-jet+7 cross-section results [^8| disfavour the 
virtuality ordering employed in JETSET's parton shower [f21 j. This has led to 



some development of the Monte Carlo shower algorithms |p7| . A more critical 
test will be the rate of soft, wide-angle photons, but here one must be especially 
wary of non-prompt photons arising from decays of particles produced in the 



hadronization 59 . 



2.5. Colour coherence phenomena 

Long-wavelength quanta see event structures on larger scales and so are 
sensitive to the presence of neighbouring charged partons M. This results 



in the effective radiating units being charge-anticharge dipoles |6(| and leads 
to a suppression of soft gauge quanta due to the requirement of coherence. 
This is true of QED and QCD but in the latter, because gluons carry colour 
charge, one must also allow at leading order for changes in the colour antennae 
formed by the hard partons. 

Effects due to colour coherence are expected in both the perturbative 
and non-perturbative stages of an event ||61|| . In the parton shower they can 
be simply incorporated by requiring angular ordering - - that is successive 
branchings are nested [§, this is known as the modified leading logarithm 
approximation (MLLA). As noted earlier, in JESTET angular ordering is 
imposed on an initially incoherent shower plj, in HERWIG it is built into 



the choice of evolution variable [17], it is intrinsic to the colour dipole model 
used by ARIADNE |12|] and it is not included in CO JETS jES]. Hadronization 
models also respect an event's colour structure in accordance with Local Parton 
Hadron Duality (LPHD) [p3"| , discussed below. Strings may be regarded as the 
natural limit of an evolving chain of colour dipoles whilst clusters may be 
thought of as the final, colour neutral, dipoles. In all Monte Carlo models the 
treatment of both stages is based on the large N c approximation to the colour 
flow |M| which gives very good agreement with known full analytic results |65j . 
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2.5.1. The string effect The classic direct test of inter-jet colour coherence 
is the 'string effect' - - a comparison of the particle flows between the jets 
in qqg events p6[. As the historic name suggests, the relative depletion 



between the qq and qg jets was expected to have a non-perturbative origin [6j 
later, assuming LPHD, a perturbative explanation was found pD |. To 
describe the Z° data [[43], |68|, it is necessary to employ a coherent shower 
with an approximately equal contribution coming from the non-perturbative 
hadronization model (but see A second classic measurement |7D| , 

proposed in f7]J, is the relative depletion between the qq pair in qqg and 



qq7 events |69], [72|, [73[| . The azimuthal angular dependence of particle flow 



in qqg events has also been studied [|74[]. In order to avoid having to find 



jets, the energy-multiplicity-multiplicity correlation can be used to measure 
the inter-jet coherence ||43||, though in this case the pQCD calculation suffers 



large corrections [75]. (A related quantity, the asymmetry in the particle- 



particle correlation, has even been suggested J76J as a way to measure intra-jet 
coherence effects |43|, [FT].) Again all of the above studies at the Z° confirm that 
Monte Carlo models incorporating a coherent shower give the best descriptions 
of the data. Unfortunately a test for the 'negative qq dipole', expected in 
qqg events when the large N c approximation is not used, does not appear 
practical |6( 



2.5.2. Fully inclusive momentum spectra The intra-jet effect of colour 
coherence is to limit the production of soft particles in the parton shower, 
leading to a 'hump-backed plateau' shape for the particle spectrum in the 
variable £ = ln(l/x) || 781 , where x is the ratio of the hadron momentum to 



the beam energy. By contrast, in the incoherent tube model, based on (||), 
there is no gluon suppression and no broad peak forms. The actual calculation 
is of the parton (essentially gluon) spectrum in a jet at the cut-off scale: 
cr _1 dcr/d£ = /(£; Q, Q , A), with Q ^> Qo^Aq C d, which it is then argued also 
applies to the hadron spectrum. The result is a rather unwieldy expression so 
that often LPHD is invoked, in the rather technical sense of taking the limit 
Qo — > A |H, to obtain the simpler, more fully evolved, limiting spectrum. 
Finally in the peak region a distorted (downward skewed, platykurtic) Gaussian 
function Jn| is applicable. The Z° data |5j|, |80], |83]. Q are qualitatively 

well described; quantitatively the limiting spectrum is a little too narrow in 
the peak region, where the distorted Gaussian fits better, although not in the 
tails. The coherent Monte Carlo event generators |12, 17, 21] give the best fits 
(however see ]|85jj). 

A particularly interesting measurement on the £ spectrum is that of the 
peak position, £*. This occurs at low momentum where the occurence of 
successive parton branchings should lessen any dependence on the primary 
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quark flavour. (Measurements of the width, skewness, etc are statistics 
limited and sensitive to primary quark flavour effects in the distribution's 
high momentum tail). In pQCD the Q dependence of £* is predicted to be 
linear in ln(Q/Qo) whie the width grows as ln 3/2 (Q/Q ) Specifically 
= a + n\n(Q/Q ) (i.e. x* = C(Qo/Q) n ) with n = 1 in the incoherent double 
logarithm approximation and n = 1/2 in the coherent MLLA, corresponding 
to a harder spectrum. The Z° data when combined with lower energy results 
SE|| ) clearly prefer n = 1/2 especially when higher-order corrections are 



; sec 



included [£9[ [87]. Similar conclusions have been reached from Breit frame 
analyses of deep inelastic scattering ep events at HERA |[8] . 



A closely related measurement is that of two-particle momentum 
correlations f89fl . Qualitative agreement is seen with the expectations based 
on the pQCD distribution of gluon pairs in a shower [pOfl . However to 
obtain quantitative agreement large corrections to the prediction should be 
anticipated |JI . 



Finally, scaling violations are seen when comparing the fragmentation 
function, <T -1 dcr/dx, measured at the Z° to those from lower energies |83|, |8~4fl , 
after making allowance for the varying primary quark favour mix (see figure [l]). 
As in deep inelastic scattering these variations are controlled [RU1 by Altarelli- 



Parisi equations (0) which allows a s to be determined |92 



2. 5. 3. Charged particle multiplicity The integral of the momentum distribu- 
tion of all particles gives the event multiplicity. Soft particles make a significant 
contribution to this total so that again colour coherence is important. The yfs 
dependence of the first few moments of the distribution of event multiplicities 
is calculable in MLLA [93||. The effect of coherence is to slow the growth of 



the mean multiplicity with ^fs |M as compared to the LLA result | 95fl . The 
MLLA prediction for the mean works well when higher-order corrections are 
included |8l], [87], |96], |97|, (5^, |99| , as do several more phenomenological functions 
of \fs |100|| , including a simple statistical phase space model ||101]| . However 
this is not the case for the incoherent COJETS model. The width (or equiva- 
lently the second binomial moment) is larger than the data |97|, |102|| although 
further relatively large higher-order corrections might be anticipated from the 
ratio of leading to next-to-leading order predictions [p5| . At high energy the 
ratios of these moments are expected to become a/s independent, a feature 
already present in the data. This would imply KNO scaling [ |103|| which is 
indeed seen to work well. The ratios of moments as given by QCD behave ap- 
proximately like those of the negative binomial distribution (NBD) ||104|| (with 
parameter 1/k « 0.4 — 0.9^/aJ), though empirically the discretized log normal 
distribution (LND) ||105|| is a better fit to the multiplicity distribution. 

Multiplicity distributions have also been studied in restricted rapidity 
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intervals |98], |106|| where the influence of global conservation laws may be 



reduced. The shape of the distribution becomes narrower for small rapidity (y) 
windows, and some structure, possibly attributable to multi-jet events |107|| , 



arises at intermediate y ranges. These features are described well by JETSET, 
but less well by HERWIG or the simple NBD and LND. 

In addition to studies of the natural flavour mix, light and heavy 
quark initiated events have been investigated separately | |108| , |109| |. A large 



quark mass reduces phase space and shields the the collinear singularity 
associated with forward gluon emission, causing a relative suppression of 



forward hadrons [0. In practice observation of the reduction in associated 



multiplicity due to this 'dead cone' jllOfl is hampered by the presence of the 



heavy hadron's decay tracks (5.5±0.1 per b-hadron ||109|| ) which actually ensure 
a higher multiplicity than in u,d,s events. However the difference in multiplicity 
is predicted ||111|| and seen to be a/s independent [p. 081 , |109| , |112|| . (A variation 



of 0.4 tracks is expected from the change in the heavy quark fraction between 
12 (T(4S)) and 91 (Z°) GeV fllQ9| ; this is below the level of the experimental 
errors) . 

2. 6. Local parton-hadron duality 

Implicit in the above calculations of event features, including those by 
Monte Carlo methods, is the assumption that pQCD provides the dominant 
contribution. This is especially true for infrared (soft gluon) sensitive 
quantities. In other words hadronization and resonance decays causes little 
disruption of the features already established by the cascade. In the absence of 
a well developed theory of hadronization this is in fact a minimum requirement 
for pQCD calculations to be worthwhile. In terms of the earlier space-time 



picture, described in subsection 1.2, this idea of 'soft confinement' appears 



very natural and it is indeed in many instances supported by data ||113|| . 
Since in this picture it is the relatively soft gluers, following in the wake 
of the hard quarks and gluons, which cause hadronization to occur, then 
except for a 'collective' action by the gluers, large momentum transfers during 
hadronization are precluded. 

An important observation which allows this picture to be taken further is 
that after a coherent shower, hadronization should occur locally. This is seen 
by considering two separating perturbative partons, with opening angle 6, and 
two gluers each emitted at 9 e to their parent partons H . At the time the two 
gluers simultaneously hadronize, their transverse and longitudinal separations 
are given by: 

6 6 
d± « R-^ > R dii » R9.— > R9 ~ R (10) 

e e e e 
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Here the inequalities follow from the strong angular ordering condition 8 € <C 
['52]. This implies that the hadronizing gluers form only a rather low density 
system in configuration space. The large inter- gluer separations effectively limit 
the influence which one hadronizing parton system can have on a neighbour. 
This locality hypothesis also receives strong support from the pre-confinement 
property of pQCD ||114|| (see section |3.3| ). 

In this scenario the local, and essentially independent, nature of each 
parton's hadronization leaves little scope for long-range effects. Therefore 
hadronization should not significantly alter an event's angular structure or 
its energy and multiplicity distributions or their correlations, but only provide 
a 'correction'. It is beyond this approach to predict actual production rates 
of specific hadrons; however it would support the notion that these rates are 
essentially constants, dependent on at most a few nearest neighbour partons. 
To a certain extent present models do respect this idea and it is the quantum 
numbers of a few nearby partons which determine the properties of a produced 
hadron. This in turn gives rise to an approximate local conservation of flavour, 
baryon number etc. 

This idea has become known as local parton-hadron duality and represents 
perhaps the simplest working hypothesis for the effects of hadronization. 
Experimental results on 'hard event properties' clearly indicate that the 
concept contains a grain of truth but also that it is manifestly untenable at 
the level of 'one parton one hadron'. Its failure should be interpreted as a need 
for a less trivial, non-perturbative model of hadronization. A more specific, 
technical definition of LPHD in terms of the limiting momentum spectrum 
(Qo — > A) IH was encountered in section 2.5.2 and will be discussed further 
in section 



3. Models for hadronization 

So far it has been argued that the parton to hadron transition should occur 
locally in space-time, only involving a few neighbouring partons, and that 
local parton-hadron duality provides a reasonable summary of its effects on 
perturbative event structures. However this does not mean that hadronization 
is simple in detail. In practice many factors influence the production of a 
specific hadron and a realistic description is presently feasible only within the 
context of stochastic, non-perturbative models. 

3. 1 . Independent fragmentation 

Independent fragmentation is perhaps the earliest framework for hadroniza- 
tion |p.l5|| , later becoming synonymous with the work of Field and Feyn- 
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man As the name suggests, the hadronization of each individual parton 
is treated in isolation as a sequence of universal, iterative qi — >q2+h branch- 
ings based on the excitation of (di) quark pairs. Unfortunately the scheme 
has no strong theoretical underpinning and is rather arbitrary in its details, 

|TT§ for 



leading to many variants. It is used by COJETS |l5fl (and ISAJET 
hadron-hadron collisions) and is a available as a set of non-default options in 
JETSET PT[. 



qi :P + 




h = (qiq 2 ) : xP + ,p T , m 2 T /xP+ 



(l-x)P + ,-p T ,p 2 T /(l-x)P + 

Figure 4. The 'unit cell' that is iterated in the independent fragmentation 
scheme, showing the light-cone momentum fractions: p+,pt,P- where 
p± = (E ±p z ). Note that the parent quark acquires a mass (m^ + p^)/x + 
p 2 T /(l-x). 

Figure f| illustrates the 'unit cell' used in the iterative implementation 
of the scheme. The flavours of the (di)quarks generated are selected in 
fixed ratios. Empirically it is found necessary to suppress both strange 
quarks and diquarks, u:d:s = l:l:r s with r s « 1/3 and q:qq' = l:r qq with 
r qq « 1/9. Further rules are required to choose between the various (low 
lying) hadron states of a given flavour. The light-cone momentum fraction, 
x = p^/p^ = (E + p z )h/(E + p z )q 1 , of the produced hadron, h = qiq2, is 
given by a longitudinal fragmentation function, such as 

f(x) = l-a + a(l + 6)(l-x) fe (11) 

In COJETS a dependence on the mass of an emitted light hadron is also 
built in. Special forms are employed in the case of heavy quarks where harder 
momentum spectra are expected (see section |5.4.2|) . The transverse momentum 
is chosen from a Gaussian distribution exp(— p\.j<y 2 ) possibly with a width that 
narrows as x — > and X ► 1 5 clS expected from phase space considerations. 
The iteration continues until a backward moving hadron would be produced 
(p z < 0) or, in the case of COJETS the sum of jet masses violates an 



available energy bound |[L17|| . Diquarks are treated just as quarks, whilst gluons 
are first split into a light qq pair (the momentum is shared equally in COJETS) 
and the above algorithm used but with retuned parameters. 

Typically a large cut-off value, Qo ~ 3 GeV, is used, resulting in only 
a few final-state partons. Since these partons are treated in strict isolation, 
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essentially ad hoc remedies must be used to ensure global conservation 
of quantum numbers. Conserving four-momentum after the partons have 
acquired masses proves particularly troublesome because event shape variables, 
and hence a s determinations, are sensitive to the nature of the chosen 
solution A fully Lorentz-invariant scheme has been proposed in ||119|| but 

it is hard to implement in general and only JETSET contains a simplified 
version, as an option. 

In addition to the independent hadronization of the final-state partons, 
independent fragmentation models also naturally employ incoherent parton 
showers, the combination of the two features offering, by today's standards, 
a mediocre description of the exacting Z° data. However when applied to 
the hadronization of a back-to-back qq pair there is little practical difference 
between independent and string-based models. (It is only after gluon jets occur 
that differences become apparent). One can therefore speculate that a possible 
way to improve the situation is, after the forced gluon splitting, to apply the 
model to each neighbouring, colour singlet qq' pair in its own CoM frame just 
as to a pair of back-to-back jets. 



3.2. String models 



When a pair of oppositely coloured quarks move apart it is thought that the 
self-interacting colour field between them collapses into a long, narrow flux 
tube/vortex line, called a string. Neglecting a short-range 'Coulomb' term the 
energy of this system appears to grow linearly with the separation. That is, 
the string has a uniform (rest) energy density or constant tension, estimated to 
be k pa 1 Gev/fm. This is equivalent to a linear confining potential as expected 
from Regge phenomenology (k = l/2ita') | 120 |, bag model calculations ||121|| , 
lattice studies |p.22j| and quarkonium spectroscopy ||123| . This picture of a 



collapsed field is analogous to a chain of magnets ||124j| and the behaviour of 
magnetic fields in (type I/II) superconductors ||125 . 



The transverse size of a string, (r 2 ) = it /(2k), is small compared to a typical 
length. Therefore it is reasonable to try to model the string dynamics on those 
of a 'massless', relativistic string possessing no transverse excitations. The 
classical equations of motion, derived from a covariant area action, imply that 
in the CoM frame the two end quarks simply oscillate backwards and forwards 



along a line in what are known as yo-yo modes ||126|| , seen as diamonds in 
figure 0. The equations of motion also admit solutions that include localized 
energy-momentum carrying 'kinks' which have successfully been identified with 
hard gluons [127|. At the end of a perturbative shower, string segments 
develop between neighbouring partons, the ends terminating on quarks. The 
full three-dimensional treatment of such a string system is rather complex, 
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being characterized by the appearance of new string segments as intermediate 
gluons lose all their momentum to the system. However, robust, covariant 
evolution algorithms are available ||128| , 129|| , 



and whilst ambiguities exist 
they are largely confined to the matching at gluon 'corners'. The classical 
equations of motion are also more complex, becoming non-linear when the end 
quarks are massive |130|| . It is worth noting at this point that the quantized 
theory of the idealized string has spawned a rich and still growing subject 
of its own ||131||, which may still prove of relevance to the hadronization and 



confinement problems [|132 



An alternative scenario for a string description of gluons arises when the 
possibility of an octet colour flux tube is admitted. For example, in a qqg 
event the quarks may be attached to triplet strings and the gluon to an octet 
string, all three of which join at a central node. A bag model calcul ation ||121 
suggests the ratio of octet to singlet string tensions is r = k 8 /k 3 = JCa/Cf 



If 

however r > 2 (as suggested by lattice calculations) then it becomes favourable 
for the octet string to split into two triplet strings and the above picture 
is recovered. The Montevay independent fragmentation model ||119|| requires 
choosing a frame in which the central junction is at rest. 

The above describes the motion of an idealized classical string due to the 
exchange of energy between the end quarks and the string; in reality a second 
process is also believed to contribute. Quantum mechanical effects allow the 
the creation of qq or qq'qq' pairs in the colour field of a stretched string, causing 
it to break in two a la the snapping of a magnet. This picture is the basis for 
the Lund group's familiar string hadronization model [jrj|, |128| , |133| |. However 



the Lund version is only one of several possible [|134| , |135|, [136 , as illustrated 
in figure |5|. 



String area law 



clusters 



breakup products 



mass regulator 



Artru-Mennessier 



quarks 



Cal-Tech II 



hadrons 



LSFF 



flavour rules 



Lund 



hadrons 



UCLA 



Figure 5. A family tree for string models. 



The starting point for all string-inspired hadronization models is Wilson's 
exponential area decay law, dV/dA = P exp(-P A) fl22j . This describes the 
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probability of a string break occurring due to the creation of a qq or qq'qq' 
pair in the colour field at a point containing the space-time area A within its 
backward light cone. As strings are believed to be uniform along their length, 
the probability of pair creation, P , is a constant per unit area. Now the 
momentum of a string fragment is proportional to its spatial extent (E = kAx 
and pz = kA€) such that the string fragment mass is given by m. 2 tring = 2k 2 A, 
so that the decay law can be reformulated as d"P/dm 2 = 6exp(— bm 2 ) with 
b = P /(2k 2 ). Since (A) = Pq 1 this implies that on average the string break-up 
points, and hence the hadron formation points, lie scattered about a hyperbola: 
t 2 = t 2 — x 2 = 4/Po ~ 2(m 2 tring )/fi: 2 . This in turn implies that the slowest 
moving fragments form first near the centre of the string (this is true in any 
frame as the break-up points are space-like separated) — an inside-out pattern 
is assured 0. 

Before constructing an actual model one must decide what the string 



fragments are to be identified with: either continuous mass fragments [134, 



137 1 - clusters — which then decay into hadrons, or actual discrete- mass 



hadrons ||138|| . The first choice leads to the original Artru-Mennessier [ |134|| 



or CalTech-II [|135|| schemes and the latter to the Lund |L4] and UCLA [|136|| 
schemes. 

3.2.1. The Artru-Mennessier /CalTech-II schemes Repeatedly applying the 
area decay law alone will result in an infinite sequence of ever smaller string 
fragments ||134| , |139|| . One way to see this is to note that the area law 



is equivalent to a joint distribution in light-cone momentum fraction x and 
transverse mass m^, which reduces to a divergent m\ distribution: 

d 2 P b ( ,m 2 i\ dV . 2 m 2 ± ->o 2 

bEi(bm\) ~ ln(6mj_) (12) 




dxdm 2 x \ x J dm^ 



(Ei is the first exponential integral function.) Since a physical interpretation 
of string fragments with very low mass is implausible, in practice a regulator 
is required. 



In the more fully developed CalTech-II model ||135|| this regulator is 
supplied by introducing a probability to allow a given string fragment to 
undergo any further splitting: 

(^string - ^o) 2 " 



V (further break) = 0(m string — m ) 



1 - exp -P - 



2k 



2 



(13) 



where Pq and m are free parameters, the latter related to the threshold mass 
for a string's decay to two hadrons. The function is the Heaviside step 
function: Q(x) = for x < and Q(x) = 1 for x > 1. If allowed, a break-up 
point is selected according to the area law with the qq or qq'qq' flavour chosen 
according to fixed probabilities from those kinematically allowed; at this point 
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no transverse momentum is introduced. Occasionally a final string fragment 
is below a second cut-off and this is replaced by a single hadron. Otherwise 
the fragments undergo a comparatively complex sequence of cluster decays in 
which phase space determines the produced flavours and momenta. 

Whilst the CalTech-II model has attempted to combine the desirable 
features of both string and cluster models, its success in confronting data has 
been at best mixed. Consequently the Monte Carlo program has fallen out of 
favour and is no longer being actively developed. 



3.2.2. The Lund scheme The alternative to a continuous mass spectrum for 
the string fragments is a discrete spectrum, the allowed values of which are 
identified with the masses of known hadrons. This approach is followed by 
the Lund ||138|| and UCLA [|136|| models. Requiring a string to fragment into 
hadrons of given (transverse) mass leaves only the choice of the hadrons' light- 
cone momentum fractions, x, free. These x values may be iteratively chosen 
according to several possible distributions and still remain true to the 
area decay law. However a set of plausible assumptions greatly reduces the 
number of allowed fragmentation functions [[1401 , |141|| . The assumptions are: 



the equations of motion are those of a classical, relativistic, constant tension 
string with no transverse excitations; a statistical left-right symmetry; a central 
rapidity plateau; and negligible end effects. The resulting x-distributions are 
known as the Lund symmetric fragmentation functions (LSSF). 
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(14) 



The coefficient a a relates to the parent quark flavour and ap to that of the 
quark or diquark produced in the colour field: in practice only diquarks are 
allowed a different value of a [^TJ. Taking every a a = a the LSFF simplifies 
(compare to (0) where a = 0) and is equivalent ||140|| to a flat rapidity, 
y = l/2\n[(E + P Z )/(E — P z )], distribution and the Wilson area law modified 
by a perimeter (Coulombic) term. 

An important issue now is to prescribe how the actual hadrons are chosen, 
and here again the models diverge. The Lund approach [|142j| is based on 
an attempt to model, principally through flavour and spin selection rules, 
the supposed quark dynamics in the strong colour field that is a string. 
The idea is a development of the concept of fermion pair production in a 
strong electromagnetic field [ 143 , 144]. To supply the energy for a qq pair, 



each of transverse mass m 2 < 



m; 



+ p±, it is necessary to consume a finite 



length of string (2ra_i_//c). The quarks have equal and opposite p± since 
no transverse string excitations are permitted in the model. If the qq pair 
are produced locally at a point then they must tunnel out to this classically 
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required separation. Using the WKB approximation for the matching of the 
quark wavefunctions at the classical turning points suggests a suppression 
factor P3|,|I31: 

exp ^(rriq = exp ( — ~ m ql x ex P ( — ~P±1 (15) 

This would be a crude approximation to the known full QED expression for 
the production of a single qq pair, I] n n 2 exp(— mrm 2 ± /n). Unfortunately it is 
not known what to use for the quark masses in (|T5|) (see reference [[[]] for some 
discussion of the range spanned by current and constituent quark masses) and 
so only qualitative conclusions can be drawn. Among these are: 

• The transverse momentum suppression is the same for all quark flavours 

• Charm and bottom quarks will not be produced from the string 

• Since m s > m u « m d (SU(3)p is broken), strange quarks will be suppressed 
relative to up and down quarks: 

= 7s < 1 (16) 

• Production of diquarks will be suppressed: 

^M<1 (17) 
7>(u) < P(uur 7s P(uu) 7s 1 } 

the power 4 arising in the last term because of the quadratic dependence 

of the tunnelling probability on the diquark mass. 

Also the Gaussian ^-distribution in (ff3j) is best regarded as a first 
approximation, leaving open the possibility of long tails originating from 
unresolved gluon emission [21, 145|| . This may account for the fact that the 
prediction for the width of the hadron p^-distribution, Jk/tt = 0.25 GeV, 
proves too narrow compared with measurements of ~ 0.40 GeV f|3], fPfl . 

As yet the model does not supply as much guidance on how to account 
for any possible spin dynamics. Two factors influence the relative production 
ratios of same flavour mesons [[14], |146|| . First are the (2 J + 1) spin counting 
factors. Second is the need for the quark produced in the string to match onto 
the wavefunction of the produced meson. At the classical boundary to the 




tunnelling region the meson wavefunction is expected to behave as 1/ ■Jm±_ [14 
The result for the vector to pseudoscalar meson production ratio is therefore: 

^Kp)J (18) 

As (my — mp)/(my + rnp) is 0.69 for ir/p and 0.004 for B/B*, the primary 
p/n ratio is taken to be 1 whilst the B*/B ratio is fixed at 3. Radially and 
orbitally excited mesons are expected to be suppressed pTl, although L = 1 
mesons may be included at the expense of a new parameter for each of the four 
additional families of states. 
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The fact that the transverse mass is involved in the suppression factor 
proves unimportant except for pions, since (p 2 ± ) ~ (0.3 GeV) 2 ^m^.. Hence it 
is anticipated that (directly produced) pions will be enhanced at low p± and 
have a tighter p± distribution than, say, p mesons or kaons. Also, because of 



their small mass, neighbouring-pion correlations are anticipated [IT?]. 

To summarize, the Lund prescription for an iterative string fragmentation 
scheme is: first choose a quark flavour; then choose a produced hadron species; 
select a quark transverse momentum; and finally select an x value. The 
normalization, N, of the LSFF in ( |14|) for this case is thus: 



AT" 1 = J dx (1 x X)<1 exp [-b^pj = F(m 2 ± ) (19) 

The resulting model requires the specification of a relatively large number of 
free parameters . Much of this can be traced back to the unknown properties 
and dynamics of (di) quarks. The UCLA model [|136|1 attempts to finesse these 
problems by formulating an iterative scheme only in terms of known hadron 
properties, thereby trying to avoid the issue of quark production in a string 

3.2.3. The UCLA scheme At the heart of the difference between the Lund 
and UCLA models is a reinterpretation of the LSFF [ 148|| . In the UCLA 



model [|136| the LSFF is used to choose both x and the species of produced 



hadron. This means that iV (|l9|) becomes an absolute normalization: 



iv- 1 

h 



£(CG) 2 JdplF(ml+pl) (20) 



The sum runs over all hadrons containing the parent quark, with CG the 
appropriate Clebsch-Gordon coefficients for the hadron wavefunction. Since 
N is now a common constant the hadron mass dependence appearing in the 
exponential term in ( |1~4"D immediately implies a suppression of heavy hadrons 



and a stiffening of their fragmentation function ||148[ |. 

An instructive way to view the difference between the Lund and UCLA 
approaches is to consider the complete weight for the production of a set of iV 
hadrons {hj} in a string of mass s stretched between two quarks qo and q^, as 
illustrated in figure |5|. 

The 'master' equation for this process's weight is given by: 

N N-l 

dW^(s) = *W(P Q - J2p hi )eM~bA} x J] exp{-«(pi) q J 

i=l i=l 

N 

x niVhJCG(q i _ 1 ,q i ;h i )] 2 exp{-/ t '(pi)h 1 }d 4 p h ^(p^ -m 2 J (21) 

The first two terms impose overall four-momentum conservation and the area 
decay law. The second set of terms is associated with the production of 
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Figure 6. A schematic diagram of a fragmenting string: the Si control 
the production of quark flavours; the "knitting" factors Ni control which 
hadrons form from the quarks. 

the N — 1 intervening qq pairs: if> qi is a possible quark flavour suppression 
factor and k controls the quark's Gaussian p_i_-distribution. The third set 
of terms is associated with the N hadrons {hj} formed out of the quark- 
antiquark pairs {(qi-iqi)}: CG is the Clebsch- Gordon coefficient for the 
hadron's SU(6) wavefunction, is an additional "Knitting factor" and k' 
controls the hadron's assumed Gaussian ^-distribution. The approaches to 
these various terms in the UCLA and Lund models are summarized in table 1. 



Table 1. The two sets of factors appearing in (21) and figure 6 in the UCLA 
and Lund string models indicating the typical values assigned. 



Factor 


UCLA [ 


136| 


Lund HD 


Controls 


S q 

K 
K 1 

N h 


1 


(n-2)6/(n+l) 
N (const.) 


7 S , (ud)i/(ud) , . 

"tt/k" 



F(mi) -1 {VyP, . 


quark flavour suppression 
quark p\ suppression 
hadron p\ suppression 
. .} hadron "knitting" factor 



A number of points in the iterative implementation of the UCLA scheme 
are noteworthy. First, the original model was designed to look ahead to the 
next iterate. For example, if the first hadron leaves behind a u or an s quark 
then the next hadron is most likely to be a pion or a kaon, and the latter choice 
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is (doubly) suppressed because of the higher masses of strange hadrons. This 
was an attempt to mimic the quantum mechanical projection of the string onto 
a set of hadrons {h}. In fact, using (21) and the UCLA ansatz for an event's 
overall weight it is possible to derive an iterative scheme that automatically 
generates chains of hadrons according to this overall distribution. That is, the 
projection from partons to hadrons is automatically taken into account. The 
required fragmentation function is remarkably similar to the original LSFF: 

exp I _ b ,<+PL/(n-l))\ (22) 




X \ X iJ J \ X 



Compared to QHD a finite mass correction term appears together with a 
term, p\/{n — 1), coming from local transverse momentum compensation, 
(n pa 2). In the case of hadrons containing a heavy quark, x is replaced 
by x e $ > x which softens the momentum spectrum and favours heavier 



states ||136|| . Second, 'multiple' so-called popcorn (discussed below) baryon 
production BB, BMB, BMMB ... is included (here B represents a baryon and 
M a meson), but due to its slow 'convergence' an ad hoc suppression needs to be 
introduced: exp(— rj m M ) with rj free (~ 10 GeV -1 ). Third, the mechanism 
for the local conservation of hadron transverse momentum proves awkward due 
to an ambiguity between the quark and hadron levels. Finally, but perhaps 
most significantly, the model only contains four (+two) free parameters: a, b, n 
and rj (+A and Qq). 



3.2.4- Consequences of the string's space-time structure A further interest- 
ing aspect of string models is that inferences can be drawn from their associated 
space-time picture. These include predictions on spin correlations and effects 
due to quantum statistics. 

If a qq pair are produced with some (equal and opposite) transverse 
momentum, p±, with respect to the string, then because they are separated by 
a distance 2m±/ n a non-zero angular momentum L = 2m±p±/K is necessarily 
introduced. Since total angular momentum J = L® S q ® Sq must be conserved 
and (L) ~ lh the qq pair typically form in a 3 Po state, particularly so at 
higher p±. This is expected to lead to p±— transverse-spin correlations, and 
spin correlations between neighbouring hadrons ||149|| . Transverse polarizations 
are only possible because the string introduces a preferred axis. (Such an axis 
is implicit in the chain-like structure found in cluster models, though it is 
not presently utilized in them). Conservation of angular momentum is also 
expected to lead to a suppression of orbitally excited hadrons. 

If two identical hadrons are produced from a string there exists an 
ambiguity in the possible rank ordering of the particles, unless they have 
exactly equal momentum. These two string drawings, illustrated in figure [7], 
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Figure 7. The two sequences of string breakings possible when two identical 
particles are present; all other particles are the same. This results in a 
difference AA between the enclosed total string areas. 



enclose different areas. Now if, as is believed, the string area law derives from 
the modulus squared of a matrix element M. = exp[L4(K+iP / '2)] then quantum 
mechanical interference will occur [ |150|| . Dependent on the Bose-Einstein or 
Fermi-Dirac nature of the particles the expected joint production probability 
becomes 

cos(ftAv4) 



(Pl2 + P 



21, 



1 ± 



cosh(P AA/2) 



(23) 



which is clearly modified from the naive sum of the weights. (A quark 
transverse momentum correction, (it/2K)Ap\, is not shown in (|23|).) As the 
momentum difference squared between the identical hadrons vanishes, AA — > 
and any enhancement or cancellation becomes maximal. The area difference, 
AA, is related to the size of the emitting volume, likened to the Hanbury-Brown 
Twiss effect in optical astronomy [151], but not related directly to the total 
string size. This approach therefore has close parallels to a more geometric 
picture based on the Fourier transform of the source distribution | 152| |. Only 
identical neutral pions can be produced side by side from a string — identical 
charged pions must have at least one intervening hadron. Therefore AA will 
be larger for charged pions, leading to smaller correlations than for neutral 
pions ||145|| . The effects of correlations on short lived resonance decays can also 
be included ||153| , |154| . 

Two basic algorithms are currently available for including a Bose-Einstein 
event 'weight' []155|| . The standard scheme available in JETSET [ [2~T| |156[ 
involves rearranging identical boson momenta so that they are distributed 
according to the pairwise correlation function. Full multiboson correlations 
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may be included, but at the cost of additional computing time. The effect of 
including Bose-Einstein correlations may be likened to adding an attractive 
inter-boson force, leading to 'lumpier' distributions. The experimental 
situation is discussed in section [7T| . 



3.2.5. Baryons and the popcorn mechanism In the Lund string model, 
baryon production poses particular problems of principle. The basic difficulty 
appears to stem from incomplete knowledge of the internal structure of a 
baryon ||157|| . Is it a quark-diquark system or a three quark system? This 
ignorance poses less of a problem for HERWIG and to some extent the UCLA 
string model because these essentially only need to know a baryon's mass and 
spin. However the Lund string works directly with the (di) quarks themselves 
and so in the absence of a guiding principle it therefore allows for two options, 
the diquark | 158 | and popcorn [ |133|| baryon production mechanisms. 




B 



M 



B 




Figure 8. A schematic of baryon production in the diquark model (left) 
and 'popcorn' model (right) leading to MBB and BMB configurations 
respectively. 



The diquark mechanism is a straightforward generalization of the quark 
meson production model and was the first to be fully developed |158| . However 
a stepwise quark model for the production of baryons was the first to be 
proposed [|144|1 and implemented [|159|| . Whilst this only evolved later [|133 



(and is continuing to evolve ||160|| ) into the popcorn mechanism it would be a 
misconception to regard it as especially contrived or unnatural. When a qq'qq' 
(or qq) pair is produced in a string's colour field with the same colour as the end 
quarks they precipitate a string break: the diquark mechanism ||158||. When a 



(virtual) qq pair is produced the possibility that they have a different colour to 
that of the end quarks allows a non-zero colour field to exist between them in 
which further real qq pairs could form, leading to the sequence BB,BMB etc: 
popcorn production |133 [. Perhaps not unsurprisingly the popcorn mechanism 
requires a (modest) number of new free parameters. 

The main practical consequence, so far, of introducing popcorn production 
appears to be a lessening of the phase space correlations between baryon- 
antibaryon pairs (see section |7.3j ). The fragmentation function for baryons 
is also expected to be softened [ |161[ ] (compared to a meson with the same 
transverse mass) due to popcorn production: in particular BMB sequences 
cause a suppression of leading baryons ||160|| . (This would suggest using 
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Figure 9. The mass spectrum for colour singlet and random qiq 2 clusters 
in a u-quark jet of four initial energies. In the random sample qq pairs from 
gluon splittings are excluded. 



> a q in the LSFF (|M).) The actual level of popcorn production required 



to describe data is still the subject of debate but it may be related to the a 
parameter of the LSFF |[L61|| . Interestingly a search for the expected chain 
like structures, such as correlated AK + p systems, has failed to see any positive 
evidence ||162|. Allowing three body decays of clusters is expected to have 



similar effects to introducing popcorn production in string models. 



3.3. Cluster models 

Whilst clusters initially appeared as intermediate states in the string model 
of Artru and Mennessier ||134j| the first fully fledged cluster hadronization 
models originated at CalTech ||137| , [L63|| . (Cluster-like structures also naturally 



arose in the earlier statistical bootstrap model ||164|| and multiperipheral 
models ||165|| .) Today the scheme is best known through the HERWIG 
implementation of Webber [17, 20]. This is based on the pre-confinement 
property proved for pQCD ||114 |: at the end of the perturbative shower the mass 
and spatial distributions of colour singlet clusters of partons spanned by quark- 
antiquark pairs have a universal distribution. In practice gluons remaining 
at the cut-off scale Qo are forcibly split into light qq pairs, the Wolfram 
ansatz ||166[| , so that in the planar approximation [Ell neighbouring quark- 
antiquark pairs form colour singlets. Figure [5] illustrates the resulting cluster 
mass distributions for u-quark initiated jets at four energies, showing clearly 
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the universality of the colour singlet cluster mass distribution, in contrast to 
the distributions for random clusters. 

The mean cluster mass is of order Qq, a few GeV, and for a colour- 
coherent parton shower has a spectrum falling faster than any power ||167|| . The 
distribution is independent of the initial parton type (q or g) and virtuality: this 
is in contrast to random quark-antiquark pairs. This universality is suggestive 
of the formation of intermediate 'super-resonances' which independently 
decay into the familiar resonances. Figure [II] illustrates the stages in the 
fragmentation of a parton into hadrons via cluster hadronization. 
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Figure 10. 

hadronization. 



A schematic diagram for fragmentation via cluster 



Three masses of cluster need to be treated: 

Average Most clusters have masses, m cl , close to Q and undergo two- 
body decays p0| . A quark or light diquark pair (here the diquarks 
serve only as a mnemonic device for quantum number conservation), 
<l3 c l3/ c l3 c l4Q3Q4) is selected at random, and for a cluster of flavours qiq 2 
the putative hadrons hi =qiqs/qiq3q4 and h 2 = q2q3/q2q3<l4 are formed. 
The selection is accepted or rejected according to its phase space weight: 
(2Jhi + l)(2Jh 2 + l)p( m ci, ^h 2 )@( m ci — "ihi + m h 2 ) where J and m 
are the hadron spins and masses, and p is the CoM momentum of the 
kinematically allowed two-body decays. If accepted, the decay momenta 
are selected isotropically in the cluster rest frame. This prescription is in 
accord with the OZI rule ||168 . 



Heavy A number of clusters in the tail of the distribution are too massive for 
an isotropic decay to be plausible. The criterion used is: m c \ = m qi q 2 > rrif 
where mf = m™ ax + (m qi + m^ 2 ) n and constituent quark masses are 
used. Here the label f denotes the flavour type of the cluster and m max 
and n are variable parameters of the model. By the (repeated) device 
of introducing a light quark pair, heavy clusters are forcibly split into 
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two daughter clusters whose masses are chosen according to a power law 
spectrum and whose directions are aligned with the qi — q2 axis in the 
cluster rest frame |I7|. This appearance of a preferred, colour field, axis 



is reminiscent of a string splitting. 

Light Occasionally it is kinematically impossible for a cluster to decay into 
two hadrons. Such a cluster is decayed into the lightest hadron h =qiq2 
and four-momentum transferred to a neighbouring cluster to satisfy the 



mass-shell constraint 

In the absence of a theory of cluster decay the simplest viable approach, 
pure phase space, is adopted (Okham's razor) 0, |137|| resulting in very few 



free parameters fPT| , |135|| . Here two-body dominance of the decays is one such 



simplification ||137|| . However non-trivial matrix elements are not precluded and 
parameters controlling the relative production rates of (di) quark pairs in the 
cluster decays and also weights for the hadron representations (knitting factors) 



are made available to users ([D], |137|| . The details of heavy cluster splitting 
prove important: m max (related to the available kinetic energy) influences 
the light baryon yield, n the heavy hadron yield and both influence the 
momentum spectra [Q. The one-body cluster decay mechanism is important 
for the production of rare, heavy states such as the T or B c , and also for 
describing leading particles, to which end a number of extra parameters have 
been added JlT] . 



Despite the fact that its initial purity has been somewhat corrupted by 
the need to handle exceptional mass clusters, the model remains a simple, 
compact and predictive scheme. Since clusters are typically light, limited 
transverse momentum is automatic, hadrons with non-zero strangeness and 
baryon number are suppressed because they are heavier, and the spin ratios 
of iso-flavour hadrons follow partly from the (2J + 1) factor and partly from 
the larger masses of higher spin states. The model represents a well founded 
attempt to go as far as possible with as little as possible. 



3.4- A comparison of the main hadronization models 

All the major hadronization models, cluster, independent and string, are local, 
universal, stochastic models. They are based on a small number of recursively 
applied branchings, where at each iterative step probabilistic rules are applied 
to select flavours, spins and momenta. The main features the models are 
summarized in table 2. 

Lorentz invariance and quantum number conservation are not troublesome 
for strings or clusters but are an issue for independent hadronization models 
requiring post facto adjustments. These essentially ad hoc remedies are 
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Table 2. A comparison of the main hadronization model approaches. 



Feature Hadronization Model 





Cluster 


Independent 


String 

Lund UCLA 


Principle 


very simple 


simple 


complex less complex 


Lorentz invariant 


yes 


no 


yes 


Flavour, charge 


automatic 


ad hoc 


automatic 


etc, conservation 








Mass dep. via 


hadrons 


quarks 


quarks hadrons 


Strangeness supp. 


predicted 


free param. 


restricted params. predicted 


Baryon supp. 


predicted 


free param. 


restricted params. predicted 


J p ratios 


predicted 


free params. 


restricted params. predicted 


Limited Pt 


natural 


built in 


built in natural 


Frag. fh. 


N/A 


free 


restricted by L— R symm. 


Cut-off (Q ) dep. 


significant 


very strong 


modest 


Stability 


infrared prob. 


collincar prob. 


stable 


Limitations 


massive clusters 


requires 


light strings treated 




treated like 


large Q 


like clusters 




'strings' 







not always implemented (particularly in application to dirtier hadronic 
collisions |116| ), and, worse, physical observables are known to be sensitive 



to details of the chosen solution 118 



Only the cluster and UCLA models, with their emphasis on hadron 
properties, provide succinct basic algorithms for flavour, spin and momentum 
selection. In their initial formulations the lack of free parameters gives these 
models a laudable predictive power, which is in marked contrast to the Lund 
string and independent hadronization models. Of course, poor fits to data 
could mandate elaborations of the models which in turn may then dilute their 
predictive power. The Lund model ought to have an advantage because it is 
based on a semi-quantitative picture of an underlying dynamics but this is 
undermined by the indirect measurability of the basic (di) quark parameters. 
Thus whilst the model's parameters possess a large degree of internal coherence, 
their must be disappointment at their large number and seemingly Byzantine 
complexity. To emphasize this point, 13 inputs are needed in the basic Lund 
model |^T) to describe the 13 (assuming u<->d isospin) L = light (u,d,s) 
hadrons (the number of hadrons increases to 27 if u and d are distinguished). 
Thus the flavour and spin aspects of the Lund model show little practical 
difference to the rather similar, but assumed, rules found in an empirical 
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independent hadronization model. However a significant difference is the way 
in which a string's causal structure restricts the longitudinal fragmentation 
function to a unique family of left-right symmetric functions ||140 |. No such 
restriction applies to independent models. 

The issue of stability with respect to collinear and soft gluons arises in 
cluster and independent hadronization models. Specifically, the spectrum of 
clusters is sensitive to the emission of isolated, soft gluons. This infrared 
instability may be regarded either as a serious problem or, perhaps, as a 
warning that it is important to treat the theory correctly. Perturbative QCD 
does not like isolated colour charges, which are strongly Sudakov suppressed; 
this is responsible for the fast falling tail of the cluster mass spectrum ||167|| . 
(Recall also that many other observables, such as the mean multiplicity, are 
known to be infrared sensitive.) In independent fragmentation a similar 
problem arises when one final-state parton is replaced by two parallel partons of 
equal net energy, giving a different multiplicity. This collinear instability occurs 
essentially because the two partons are oblivious to each other's presence. 

To recap, the motivation for the rules used in these models varies from the 
Q CD-inspired, complex dynamics of strings through the minimalism of clusters 
to the simple expediency of independent fragmentation. In the subsequent 
sections we shall see how well these hadronization models compare to the 
various Z° data. However since independent hadronization makes no claim to 
be based on QCD, a fact reflected in the extreme arbitrariness of its parameters, 
we drop it from further discussion. While it is not considered a viable scheme 
for describing Z° physics it does however survive in the ISAJET Monte 
Carlo program for hadron hadron collisions. 



3.5. Colour rearrangement 



Recently interest has arisen in the possibility that soft, long wavelength, gluons 
may cause non-perturbative rearrangements in the colour structures of events 
developed during the showering stage, particularly in relation to W + W~ pair 
production |L1, |16S|| . Several phenomenological models are available based on 
the three main e + e~ event generators, and are summarized in table 3. 

The rearrangement criterion in the ARIADNE-based models is a decrease 



in the A-measure ||1 73|| which quantifies the momentum-space size of a string 
system. Only the second model [|171|| should be considered as realistic for 
typical events. In HERWIG at the end of the perturbative shower pairs of qq' 
singlet clusters may be rearranged, with fixed probability (~ 1/9), if this results 
in a reduction in the quadratic sum of their space-time sizes. In JETSET the 
rearrangement criterion is based on the space-time evolution of the strings. In 
model I the spatial overlap of 'wide' flux tubes is used to assign probabilities for 
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Table 3. The main features of colour rearrangement models. 



Feature Colour Rearrangement Model 



Basic MC 


ARIADNE 


HERWIG 


JETSET 


Author(s) 


J. Hakkinen L. Lonnblad 


B.R. Webber 


T. Sjostrand S. Todorova 




G Gustafson 




V A Khoze 


Reference (s) 


[Em 







rearrangement : 




decrease cluster 


I space-time overlap 


- criterion 


decrease A measure 


spatial size 


II string crossing 


- in shower 


no yes 


no 


no yes 


- in hadron. 


yes yes 


yes 


yes yes 


- mult, re-arr. 


yes yes 


yes 


no yes 


- inter-singlet 


yes yes 


yes 


yes yes 


- intra-singlet 


no yes 


yes 


no yes 



rearrangement, whilst in model II rearrangement occurs at the first crossing 
of the 'narrow' vortex lines. Todorova's model [Ol is an elaboration of the 



original model [|172|| , allowing for multiple reconnections, including within a 
single colour singlet system, and self-interactions; this leads to string loops - 
glueballs. 

As yet these models have not been thoroughly investigated, nor their 
consequences for Z° decays found. Note that when including colour 
rearrangement the Monte Carlos must first be retuned: for example reducing 
the A-measure also lowers the average multiplicity, which must be compensated. 
Since the physics of colour rearrangement is universal in nature, effects should 
also be anticipated in all other types of hadronic interactions, including for 
example B-hadron decays (B— > J/ipX s ) and rapidity gap events [173], where 
important constraints may be found. 



4. The colliders and experiments 

4-1. The colliders 

4-1.1. LEP The LEP machine, a 27 km circumference storage ring, was 
conceived and constructed primarily as a Z° factory. In the LEP 1 mode, 
running at and near the Z° peak, a luminosity of 2.2 x 10 31 cm _2 s _1 could 
be obtained, with beam lifetimes of up to 20 hours. Until 1992, four bunches 
of electrons and positrons crossed every 22 /xs at the experiments. During 
1993 and 1994 LEP ran in a mode allowing eight bunches, with a crossing 
time reduced to 11 /is. In 1995, so-called bunch trains were introduced, with 
four trains, each of three bunches of particles, providing a further increase in 
luminosity. The beam spot in 1995, the last year of extensive running at the 
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Z°, was 250 jum in the vertical and 5 /iin in the horizontal directions. In the 
years 1990, 1992 and 1994, LEP ran on the Z° peak to produce the largest 
possible number of events. In 1991, 1993 and 1995, the centre-of-mass energy 
was scanned across the Z° peak to allow studies of the Z° line shape and other 
tests of the electroweak theory. 

4-1-2. SLC Unlike LEP, the SLC is a single pass collider, with the electron 
and positron bunches lost after each pass. The bunch crossing frequency, 
120 Hz, is therefore very much lower than at LEP. The SLC is however capable 
of providing electron beams with large polarization. In order to increase the 
luminosity, the beams are tightly squeezed before the collisions to reduce 
their cross sectional area. This technically difficult procedure results in a 
beam spot of diameter 0.5 fim in the vertical plane and 2.3/im horizontally. 
The experiment, SLD, is able to exploit this by placing precision microvertex 
detectors only 3 cm from the interaction point. Since 1989 the SLC luminosity 
has steadily improved to around 6xl0 29 cm 2 s _1 , and electron beam polarization 
values of 80% are now routinely achieved. 

4-1.3. Experimental conditions Both LEP 1 and SLC could provide their 
host experiments with clean, low background, experimental conditions. In 
addition, the relatively low beam crossing rates allowed the experiments to 
implement bias-free triggers and to collect data with essentially no dead time. 
For studies of multihadron production at the Z° peak, initial-state radiation is 
negligible. Multihadronic Z° decay events could be triggered with efficiencies 
greater than 99% and selected offline with backgrounds (mainly due to r lepton 
pair production) of less than 1%. 

4-2. The detectors 

The four LEP detectors, ALEPH, DELPHI, L3 and OPAL, and the SLD 
detector at the SLC are typical large, multipurpose, particle physics detectors 
designed to allow measurement of events over a large part of the solid angle. 
For studies of hadronization, certain features of the detectors are particularly 
relevant. A large, active tracking volume within a strong magnetic field 
allows reconstruction of jets of charged particles and of secondary vertices from 
strange particle decays; good track momentum and direction measurement give 
accurate reconstruction of systems of two or more particles; charged particle 
identification permits the study of inclusive identified hadrons; electromagnetic 
calorimetry allows measurement of inclusive 7r° and i] mesons; and precision 
secondary vertex reconstruction and inclusive electron and muon identification 
may be used to identify charm and beauty particles. 
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The coordinate systems used by the experiments define z to be along the 
beam direction, so that the xy plane is perpendicular to the beams. Then 
r and <fi are the usual cylindrical polar coordinates, and the xy plane is also 
called the rcf> plane. The angle 9 is normally the polar angle to the z (beam) 
axis. 



4-2.1. ALEPH The ALEPH detector 175 was designed to provide high 
three-dimensional granularity with large solid angle coverage for charged 
particle tracking and for calorimetry. Its silicon vertex detector, drift chamber 
and large time projection chamber (TPC), in a 1.5 T magnetic field, give 
a resolution on momentum transverse to the beam directions of <r(l/pr) = 
0.6 x 10~ 3 (GeV/c) -1 for 45 GeV/c tracks. At low momentum, the resolution 
is dominated by multiple scattering which contributes 0.5% to <j(pt)/pt- The 
silicon vertex detector permits the measurement of track impact parameters 
with an accuracy of 25 /im+95 fim/p (with p in GeV/c) in both the r<ft and 
rz planes, allowing excellent reconstruction of charm and beauty particles. 
The TPC also measures ionization energy loss, dE/dx, giving good electron 
identification to high momenta, it /K separation at two standard deviations 
(2a) in the relativistic rise region above 2 GeV/c, and K/p separation at la 
for momenta over 5 GeV/c. The efficiency to measure K° and A particles is 
about 50% at maximum, which occurs at about 8 GeV/c momentum. The 
lead/ wire-chamber electromagnetic calorimeter has an energy resolution of 
a(E)/E = 0.18/^+0.009 and an angular resolution of (2.5/7^+0.25) mrad 
(E in GeV). Along with the dE/dx information, the calorimeter gives an 
average electron identification efficiency of 65% in hadronic jets. The efficiency 
for reconstruction of 7r° mesons peaks at about 50% at 10 GeV, falling to 
20% at 2 GeV and 10% at 30 GeV. The 7r° energy resolution is around 7%, 
independent of energy. Muons are identified using the hadron calorimeter and 
muon chambers; only muons above 3 GeV/c momentum penetrate the system, 
and these are detected with an average efficiency of 86%. 



4.2.2. DELPHI A pivotal feature of the DELPHI |7| detector is the 
particle identification capability of its ring imaging Cherenkov detectors 
(RICH). These cover both the barrel and end cap regions, and have both liquid 
and gas radiators. When combined with ionization energy loss information 
from the tracking detectors, the system gives clear identification of charged 
particles over the whole momentum range at LEP 1. The tracking detectors of 
DELPHI operate in a magnetic field of 1.23 T and consist of a silicon vertex 
detector, an inner drift chamber, a TPC and an outer detector of drift tubes. 
In addition there are two forward chambers, containing planes of drift tubes, 
to improve reconstruction at low polar angles. The system gives a momentum 
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resolution cr(p)/p of 0.7% to 1.4% in the barrel region, varying with track 
momentum. Measurement of Kj? and A in multihadron events has an efficiency 
of 30-40% over a wide momentum range. The silicon vertex detector allows 
track impact parameters to be measured with an accuracy which depends on 
momentum and polar angle, and which varies between 30 and 130 fim in r(p and 
between 40 and 200 fim in z. DELPHI'S electromagnetic calorimeters consist 
of high-density projection chambers in the barrel region and lead glass blocks 
in the endcap regions. They are preceded by 0.8/ sin 9 radiation lengths of 
material in the barrel region, and more in the endcap regions, so that resolution 
is somewhat degraded. When combined with dE/dx information, the system 
allows electron identification with efficiency and purity values both around 50% 
over a wide momentum range. Photons are identified both in the calorimeters 
and by measuring conversion e + e~ pairs in the TPC, allowing reconstruction 
of 7T° mesons. Muons above 3 GeV/c are detected using a hadron calorimeter 
and muon chambers with an efficiency between 75 and 85%. 



4-2.3. L3 The L3 detector ||177|| was designed with the primary aim of 
reconstructing electrons, muons and photons. It is therefore more limited 
than the other detectors in its capabilities for studying hadronization. L3 
has a low magnetic field of 0.5 T in a large cylindrical volume of diameter 
12 m. The central time expansion chamber measures tracks out to a radius of 
31.7 cm with a high spatial resolution in the r<j) plane. It is supplemented by 
a z-chamber to measure the polar angle of tracks. The L3 arrangement gives 
optimized momentum resolution for penetrating muons, with cr(p)/p w 2.5% 
for 45 GeV/c muons. Its bismuth germanium oxide (BGO) electromagnetic 
calorimeter is preceded by less than 10% of a radiation length in the barrel 
region, and has a spatial resolution better 2 mm above 2 GeV. The energy 
resolution is about 5% at 100 MeV and 1.4% at 45 GeV. The calorimeter 
permits electron identification, with only 0.1% probability of misidentifying 
a hadron. It is also well suited for measurements of inclusive 7r°, rj and r/ 
production. 



4-2.4- OPAL The OPAL ||1 78|| detector has a warm solenoid providing a 
magnetic field of 0.435 T. The main central tracking jet chamber lies outside 
of a silicon microvertex detector and a precision vertex drift chamber; a set 
of ^-chambers around the jet chamber give precise measurement of track 
polar angles. The combination of the tracking chambers gives a momentum 



resolution of cr{p T )/p T « y0.02 2 + (0.0015pt) 2 with in GeV/c. Efficiency 
for reconstruction of K° and A particles varies with momentum, having a 
maximum value of 30% at about 5 GeV/c. The silicon detector, orginally an 
r0 device but improved in 1993 to also measure z, gives an impact parameter 
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resolution of 15 fim in rip for high momentum tracks. The jet chamber 
measures ionization energy loss, dE/dx, of tracks in multihadronic events with 
a resolution of 3.8%, allowing excellent identification efficiency and high purity 
for electrons, pions, kaons and protons over almost the whole momentum range 
at LEP 1. OPAL's lead glass electromagnetic calorimeter, which is preceded by 
some two radiation lengths of material, has an energy resolution varying with 
energy from 1 to 5%. The muon chambers, together with the instrumented 
flux return of the hadron calorimeter, are highly efficient for identification of 
muons with momentum greater than 3 GeV/c. 



4.2.5. SLD The SLD experiment, like DELPHI, makes use of 

Cherenkov ring imaging to identify charged particles, with a detector which 
uses both liquid and gas radiators to allow coverage of a wide momentum 
range. Tracking is done by a central drift chamber within a 0.6 T magnetic field 



which gives a momentum resolution, o{j>t)IVt = yO.005 2 + (0.01/pr) 2 , with 
Pt in GeV/c. For electromagnetic calorimetry, the SLD uses a liquid argon 
device, inside the magnet coil, which gives a resolution of around 15%/y/E 
(with E in GeV). The SLD has silicon charge-coupled-device pixel detectors 
for microvertex measurements. The pixels are 22 micron square, and the setup 
covers radii from 3 to 4 cm from the interaction vertex. Resolution on track 
impact parameter in rip is in the range 11 to 20 /im, depending on track 
momentum and polar angle. Muons are identified by layers of streamer tubes 
between the slabs of iron which make up SLD's warm iron calorimeter. 



5. Measurements of inclusive single identified particles 

Inclusive single identified particles are usually studied in terms of their 
fractional energy (xe) or momentum (x p ) relative to that of the beams, with the 
fragmentation functions being reported as (1 / cruder /dx. Here is the total 
cross section for e + e~— » Z°/7* — > hadrons. Its inclusion is experimentally 
advantageous since it obviates the need to measure absolute cross sections, 
so reducing systematic errors: da/a^ is simply calculated as the number of 
particles, AiV, in bin dx relative to the total number, N to t, of Z° hadronic 
decays. In a real measurement the bin width is finite, Ax, and the differential 
cross section is taken as l/iV tot x AN/ Ax. Because fragmentation functions 
vary rapidly with has to be taken for large values of Ax in interpreting a 



measurement as a differential cross section at some particular value of x 1 180 ]. 

Total inclusive particle yields, or average multiplicities per hadronic Z° 
decay, are obtained by integrating the measured fragmentation functions and 
extrapolating into any unmeasured regions of x with the aid of one or more 
models or interpolation functions. Systematic errors are included to account 
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for uncertainties in this procedure. 

Fragmentation functions and total inclusive yields have been measured for 
an impressively large number of particle species at LEP. The string and cluster 
models, as implemented in JETSET and HERWIG, are usually confronted with 
the data. A comprehensive compilation of inclusive particle production data 
in e + e~ annihilation at all available CoM energies above the T mass, as of 
mid-1995, may be found in [^] where measured fragmentation functions are 
plotted along with curves obtained from JETSET version 7.4. A more recent 
review ||181|| contains a good summary of measurements published after |86 . 



5.1. Overall inclusive rates 

Tables 4 and 5 list of all the measured inclusive yields of mesons and 
baryons published to date. Where an experiment has reported more than 
one measurement, only the most recent is taken. For each measurement of 
a particular particle, statistical and systematic errors have been combined in 
quadrature; then the weighted mean of the available measurements has been 
calculated to give the results shown in the tables (no attempt has been made 
to take into account systematic errors correlated between experiments). Yields 
reported over a restricted x range are given separately. 

In general there is very good agreement among the measurements of the 
different experiments. In only two cases, where the measurements are listed 
separately in the tables, is there evidence of disagreement: the DELPHI and 
OPAL measurements of the A(1232) ±± are possibly inconsistent, and the 
S(1530)° yield reported by DELPHI does not agree well with the numbers 
given by ALEPH and OPAL. 



5.2. Conclusions for Monte Carlo models 

In tables 4 and 5 the measured rates are compared with the outcome of the 
three major Monte Carlo models which attempt a full simulation of particle 
production. The numbers in bold font show results which are more than three 
standard deviations from the experimental measurements. In each case, the 
default versions of the programs have been used, although for JETSET version 
7.4 various sets of alternative parameters have been suggested which improve 
the agreement with the overall rates. The most recent HERWIG version 5.9 
does not fit as well as version 5.8, but a new default set of parameters will no 
doubt follow careful comparisons with data. 

It is clear from the various JETSET tunings suggested by the four LEP 
experiments in |ITJ that there are strong correlations among the program's 
parameters, possibly such that there is no unique best set. To this extent 
JETSET may be underconstrained despite the large number of experimental 
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Table 4. Average measured charged particle and identified meson 
multiplicities in Z° decay together with the rates from Monte Carlo 
models. The letters ADLMO indicate the contributing experiments. Where 
appropriate, the rates always include both particle and antiparticlc. 



Particle Multiplicity 



HERWIG59 JETSET74 UCLA74 Comments 



Charged 


20.96 ±0.18 


20.40 


20.95 


20.88 


7T+ 


17.05 


± 0.43 


16.62 


16.95 


17.04 


7T° 


9.39 


± 0.44 


10.15 


9.59 


9.61 


1 


0.282 


± 0.022 


0.246 


0.286 


0.232 




n no 
0.93 


i n nn 
± 0.U9 


0.92 


1.00 


0.78 


p(770)° 


1.29 


±0.12 


1.12 


1.50 


1.17 


w(782) 


1.11 


±0.11 


1.05 


1.35 


1.01 


t?'(958) 


0.064 


±0.014 


0.071 


0.127 


0.061 




0.25 


±0.04 


0.143 


0.297 


0.121 


fo(980) 


0.098 


±0.016 


0.068 






</>(1020) 


0.108 


±0.005 


0.181 


0.194 


0.132 


f 2 (1270) 


0.170 


± 0.043 


0.137 






£(1525) 


0.020 


±0.008 


0.021 






K+ 


2.37 


±0.11 


2.08 


2.30 


2.24 


K° 


2.010 


±0.029 


1.87 


2.07 


2.06 


K*(892)+ 


0.714 


± 0.044 


0.524 


1.10 


0.779 


K*(892)° 


0.759 


±0.032 


0.530 


1.10 


0.760 


K*(1430)° 


0.079 


± 0.040 


0.067 








0.19 


±0.07 


0.054 







adlmo n m m n 



O [182| 



ADL m [[H |18 

a m x E > o.i 



L 1 186] 



AD [187, 18S 



AL |87j, |189| 
A (43| xe > 0.1 



L 189] 



D [188| x E > 0.06 



ADO [187, 190, 191| 



D [188| x e > 0.05 



D [192| 



do mim 



ADLO |19J |l88j, |l85j, |195 [ 
ADO 



ADO |87, |90, |191| 



D [190| 



O [191| x E < 0.3 



D+ 0.187 ±0.014 

D° 0.462 ±0.026 

D*(2010)+ 0.181 ±0.010 



0.190 0.174 0.196 

0.406 0.490 0.497 

0.151 0.242 0.227 



ADO |97] 198, |199| 



ADO |97|, |98, |199| 



ADO [197, 198, 199| 



D" 



0.131 ±0.020 



0.087 



0.129 



0.130 



O [199| 



B* 



0.28 ±0.03 



0.182 



0.260 



0.254 



D [200| 



B* 



0.118 ±0.024 



0.032 



D [201 



3 /if) 0.0054 ±0.0004 0.0018 0.0050 

^(3685) 0.0023 ±0.0005 0.00097 0.0019 
Xd 0.0086 ±0.0027 0.00088 



0.0050 ADLO |0g |0|, |20|, p05| 
0.0019 DO pO^ |20 

dl poa EO 



1.4 ± 0.7x 10" 4 < 1.0 X io~ 7 2.2 x 1.8 X io" 6 [206| S(3 lightest T) 
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Table 5. Measured baryon multiplicities in Z° decay together with the rates 
from Monte Carlo models. The letters ADLO indicate the contributing LEP 
experiments. Where appropriate, the rates always include both particle and 
antiparticle. 

Particle Multiplicity HERWIG59 JETSET74 UCLA74 Comments 



P 



0.98 ±0.09 1.41 



1.19 



1.09 DO 193 182 



A(1232)++ 0.079 ±0.015 0.278 0.189 0.139 

0.22 ±0.06 0.278 0.189 0.139 



D [207 | 



O [208 | 



A 0.373 ±0.007 0.605 

A(1520) 0.0213 ± 0.0028 



0.385 



0.332 ADLO [194, 20E, 185, 21C] 



O [210| 





0.092 ±0.017 0.123 


0.073 


0.061 


xr 


0.084 ±0.017 0.102 


0.068 


0.056 


£+ + £- 


0.174 ± 0.021 0.225 


0.140 


0.118 




0.074 ± 0.009 0.093 


0.073 


0.074 


£*+ + £*- 


0.0474 ± 0.0024 0.202 


0.074 


0.074 




0.0265 ± 0.0009 0.0746 


0.0271 


0.0220 


S(1530)° 


0.0072 ± 0.0007 0.0352 


0.0053 


0.0081 




0.0041 ± 0.0006 0.0352 


0.0053 


0.0081 




0.0068 ± 0.0007 0.0352 


0.0053 


0.0081 




ADO 

a m 



212, 21C| 



D [212| 



O [210| 



0.0012 ±0.0002 0.0093 0.00072 0.0011 ADO [43, 213, 210| 



A+ 



0.078 ±0.017 0.0129 



0.059 



0.026 O [199| 



measurements. DELPHI have published comprehensive sets of tuned 
parameters for various Monte Carlo models (including ARIADNE and JETSET 
with matrix elements as well as with parton showers) in which they take account 
of event shape variables as well as inclusive identified particle rates. This 
exercise is useful but probably premature. Some of the recent measurements 
differ significantly from those used in the tuning, and many have much reduced 
errors. For example the Q~~ baryon is now known to be produced at a much 
lower rate than previously measured, and it turns out the DELPHI tuned rate 
fits better with this new rate than with the one used as input; the same is 
true of the 0(1020) meson rate. All of the models considered in ||4|| describe 



the inclusive rates reasonably well, with the exception of the performance of 
HERWIG in the baryon sector. 

5.2.1. Production of L=l mesons Although the available meson measure- 
ments are predominantly of the L = pseudoscalar and vector states, the 
presence of the L = 1 mesons shows clearly their importance in the hadroniza- 
tion; this is confirmed also in the baryon sector with the observation by OPAL 
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of the A(1520). Other scalar, axial vector and tensor mesons, together with 
orbitally excited baryons, must presumably also be produced, although large 
widths and small branching ratios will make them difficult to measure. Many 
of the lighter particles are therefore decay products of other hadrons, and care 
has to be taken in any interpretation of the data using only relative light par- 
ticle rates. So far only HERWIG 5.9 includes by default the P and some D 
wave meson states, and as the tables show their inclusion does not mar the 
agreement with data in the u,d sector. However HERWIG does poorly with 
mesons containing s quarks as well as with baryons. 



The production of the light tensor mesons is discussed by DELPHI in |fL8i 



where a comparison is made of relative rates of tensor to corresponding vector 
mesons. While the production ratio i^/ p° is similar to that for f' 2 /<j), at about 
20%, there is evidence for a lower K 2 /K* ratio, in agreement with results from 
hadroproduction experiments. This suggests an extra suppression of strange 
tensor mesons. However OPAL [|191|| measures a larger rate for and so 
the picture is not yet clear. Some evidence has also been reported [|188 



for a rise in the ratio i 2 / 'p° with meson momentum; in other words the 
fragmentation function of the tensor meson may be harder than that of the 
vector meson, as generally expected of heavier hadrons ||214|| . However, these 
are difficult resonances to measure; they have large widths, large combinatorial 
backgrounds and uncertainties in the resonance line shapes. As usual, more 
results would help. 

Only one measurement of a scalar meson, the fo(980), has been 
reported [ |188| ], with a ratio fo/p° of 0.14 ± 0.03. Thus if the f is indeed a 
conventional ++ meson (see for a mini-review of the scalars) then the 
scalar and tensor mesons are produced with similar rates. But again the 
measurements are difficult, and other studies ||187|| of the inclusive 7r + 7r~ mass 
spectrum with higher statistics have failed to report a measurement of the 
fo(980) because of systematic uncertainties. One should therefore be wary of 
too much interpretation of one measurement. 

The f (980) and a (980) have aroused interest ||215|| as potential probes 
of the Gribov confinement scenario ||216| . In this theory the QCD vacuum is 
likened to the intense QED fields expected around super charged ions, Z > 180 
(or > 137 for point-like charges) ||217 |. This results in the production of 
spatially compact 'novel vacuum scalars', identified with the fo(980), which 
are expected to be produced in relative isolation. Particular signatures include 
enhanced production at central rapidities (with respect to the thrust axis) and 
in low multiplicity events ||215|| . 

From a semiclassical point of view, the orbital angular momentum I of a 
qq pair from string fragmentation is given by < pr > xd where < pr > is the 
mean quark momentum transverse to the string and d is the size of the resulting 
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hadron [ I4"| . Typical values give I ~ 0.05h so that the rate for mesons with 
non-zero orbital angular momentum is expected to be less than about 10% of 
the rate for the corresonding L=0 mesons. The large average quark-antiquark 
separation in radially excited states also works to prevent production of these 
mesons in string fragmentation. Although the L = 1 mesons can be simulated 
in JETSET at the correct rates by the adjustment of appropriate parameters, 
the large experimental rates are nevertheless a problem for the basic model 
assumptions about production of hadrons from strings. 

5.2.2. Strangeness suppression Strangeness suppression is immediately 
apparent from both the meson and baryon measurements given in the tables. 
There is a large number of ways to determine, from the data, values for 7 S , the 
quark-level strangeness suppression assumed in the string model. Results are 



tabulated for example in |TTJ where all measurements agree on a value of 0.3, a 
result which accords with reasonable values for the strange quark mass and 
the string energy density k. Since the various methods in [[Hp use both light and 



heavy quark states, this consistency suggests that the suppression occurs at the 
quark level, in agreement with the string-model assumption. On the other hand 
the UCLA model also reproduces reasonably well the strange particle rates, and 
previous versions of HERWIG have been tuned to do so. And hadronization 
studies in ep collisions at HERA ||218|| give a lower value of 7 S ~ 0.2, in 



apparent disagreement with the Z° decay measurements. Therefore it is fair 
to say that the data are not yet conclusive. A direct method to measure the 
strangeness suppression in e + e~— >■ Z°/7* — > hadrons has been proposed 



which makes use of the electroweak forward-backward asymmetry and which 
could possibly distinguish between quark-level and hadron-level suppression. 
Recently SLD ||220|| have applied this method to their data on inclusive K*° 



and K*° production using 150k events, with the result 7 S = 0.26 ± 0.12. 

5.2.3. Relative rates of vector and pseudoscalar mesons In the absence of 
mass effects, the ratio of direct pseudoscalar to vector meson production may 
depend simply on spin statistics, in which case the value P/(P+V) would be 
expected to be 3/4. However feed-down from decays is also important and 
may obscure the interpretation of the experimental results. One approach to 
determine the underlying P/ (P+V) value is to tune the appropriate parameter 
in the JETSET model, but this can only be done within the limited knowledge 
available on production of the higher states. Alternatively, one can use the 
measurements in the b and c sectors where vector to pseudoscalar mass 
differences are much smaller and there are some hopes to measure the orbitally 
excited states. The average ratio for primary B*/(B+B*) is found to be 



0.75 ± 0.04 |200| , |221| , [222J, in excellent agreement with simple spin counting. 
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However the picture is not so clear when charm is considered: here the ratio 



D*/(D+D*) is measured [ 197 | to be 0.51 ±0.04. The difficulty here arises from 
incomplete knowledge of the production rate and decay modes of the orbitally 
excited D** states which are likely to feed down into the D and D* production. 
The question then arises as to why B** production does not similarly muddy 
the waters in the b sector; it may be that production and decay rates of the 
four different J p states of B** conspire to leave the value of P/(P+V) at 3/4. 
More measurements are needed before definite conclusions can be reached. 

5.2.4- Baryon production In the string model, baryon yields are determined 
by many parameters. The overall baryon rates relative to mesons depend on the 
relative probability to produce a diquark pair from the sea. Spin-1 diquarks 
may be suppressed relative to spin-0 diquarks. The strangeness suppression 
enters in a similar way as for meson production, but there is in addition the 
possibility of extra suppression of strangeness in a diquark. And the popcorn 
mechanism may introduce one or more mesons locally in phase space between 
a baryon and an antibaryon. 

Although JETSET does rather well, there are some discrepancies with 



measured rates. Attempts by OPAL ||210|| to tune the parameters which 
control baryon production have shown that it is not possible to reproduce 
simultaneously all of the measured rates. The suggestion then is that the 
mechanisms for baryon production in the string model, and particularly for 
the strangeness suppression, are deficient. However, as has been said, there 
is now clear experimental evidence also for orbitally excited states in baryon 
production. The rate for the J p = §" A(1520), at 0.02 per hadronic Z° decay, 
is around 5% of that for the J p = | A. And there are many similar baryon 
states which cannot be measured experimentally but which, it is fair to assume, 
must be produced in the hadronization. So since JETSET, like all of the 
other models, does not include production of orbitally excited baryons, no 
clear conclusions can yet be reached. 

Neither HERWIG nor UCLA, both of which rely only on phase space, 
mass and spin, are successful in the strange baryon sector, with the former 
consistently overestimating the rates, and the latter underestimating them. As 
presently implemented, baryon production in HERWIG does not take account 
of the appropriate SU(6) Clebsch-Gordon coefficients and this will lead to 
an overestimation of baryon production rates. As with JETSET, neither of 
these models includes the production of baryons other than the lowest lying 
L = states. In principle their inclusion would lower the predicted rates for 
primary low lying baryons, since some higher states would be produced in their 
stead. However there would be a compensatory increase in the rates due to 
feed down from decays. Since HERWIG consistently overproduces and UCLA 
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underproduces the baryons, it seems reasonable to deduce that at least one of 
them is incapable of being fixed up by this mechanism. 

The conclusion then from the baryon yields is that the string model 
of JETSET, while certainly not perfect, is in reasonable accord with the 
measurements. Recent work on the baryon production within the popcorn 
model has further improved the agreement ||160|| . Whether or not this is to 
be taken as a strong endorsement of the model is however an arguable point, 
given the large amount of freedom available to tune parameters in order to 
reproduce the observations. Both the HERWIG cluster model and the UCLA 
string model are clearly in difficulty, and it remains to be seen whether they 
can be rescued. 



5.2.5. Comparison to models of total yields A number of models have been 
proposed to treat only the overall yields of identified particles. Such models are 
necessarily of limited physical content, although they turn out to be reasonably 
successful in describing the inclusive rates. Why they do so is not at all clear. 

In the thermodynamic model ||223|| the source of particles is assumed to 
be a hadron gas in thermal and chemical equilibrium. The model has three 
parameters, a temperature, a volume and a parameter to allow for incomplete 
strange chemical equilibrium (similar to the strangeness suppression of the 
Lund string model). The model gives a good fit to the LEP data, as well 
as to lower energy e + e~ annihilation data, with a temperature of around 
170 MeV, close to Aqcd and the temperature found in the earlier statistical 
model of Hagedorn [ |164[ | . The author of the model speculates that the thermal 
equilibrium could be a feature of the quark-hadron transition, brought about by 
strong interactions. This argument was also invoked by Fermi [|101|1 to justify 
his phase space model for hadron production. However it is difficult to reconcile 
this picture with the conventional view of hadronization in e + e - annihilation 
as occuring locally in the wake of rapidly separating colour sources. 

A "striking regularity" [[224j| has been noted in the particle yields, and a 
simple formula proposed which reproduces well the observations (apart from 
the pions and possibly the fl~ baryon): N = (2J + l)/(/ m + 1) x aexp (— bM 2 ). 
Here, N is the yield for a particle of spin J and mass M, and I m is the isospin 
for baryons and a "modified" isospin for K and K* mesons and isosinglet 
pseudoscalar mesons. The introduction of the "modified" isospin appears 
rather ad hoc, although there are plausible arguments to justify it. The model 
makes no attempt to explain the yield differences between members of the same 
isomultiplet which, for example, are significant for kaons. The parameters a 
and b are fitted to the data, and the slope parameter b, at about 3.9 GeV -2 , is 
found to be the same for LEP as for lower energy measurements, implying that 
the regularity may be universal. It is unclear as to the physical origin of the 



46 



expression, and it will be interesting to see how the formula copes with future 
measurements. This observed regularity has a less successful predecessor of 
the form N = (2J+ 1) x aexp(— bM) ||225|| , which was applied to lower energy 
data. 

Another model ||226|| of the total yields is based on string fragmentation, 
and proposes a simple formula with only three parameters: a strangeness 
suppression, an effective temperature and a relative normalization factor 
between mesons and baryons. Following the string model, the rate iV of 
light meson and baryon production is taken as N = (C/Cb) x (2 J + 1) x 
(Ts)^ x ex P ( — -E'bind/T 1 )- The normalization C depends on the centre-of-mass 
energy, Cb is a relative suppression of baryons, J is the particle spin, 7^ 
gives the suppression for a hadron containing N s strange quarks, -Ebind is the 
hadron binding energy and T is the effective temperature. The model gives 
a good simultaneous fit to LEP and lower energy data, with a temperature 
of 298 ± 15 MeV and a strangeness suppression 7 S of 0.29 ± 0.02. The model 
also gives a good description of heavy flavour production. Its predictions for 
production rates of excited charm states have recently been shown to agree 
with OPAL measurements [ 227|| . 



5. 3. Rates for heavy quarkonia 

Due to their narrow widths and the availability of clean leptonic decay channels 
the principle heavy quarkonium states measured are the J /if) and T(1S, 2S,3S), 
based upon which further excitations can be reconstructed p02| , |203| , |204| , [205| , 



206| , |228|| . The production of these heavy QQ bound states is thought to be 



rather atypical of hadron production in general, especially for the T, due to 
the significant part played by perturbative physics. 

In the case of charmonium the dominant production mechanism is expected 
to be weak b hadron decays: b— >c+(cs) plus subsequent colour rearrangement. 
In the case of B Uj d mesons the J/ip branching ratios (~1.15%) have been 
previously measured at the T(4S) ||229|| , so allowing reliable predictions for 



charmonium rates at the Z°. (The presence of B s and b baryons at the Z° makes 
little difference to the inclusive b hadron branching ratio |[203|| ). Perturbative 



fragmentation contributes only at the few percent level to charmonium 
production ||230| , |231|| , as is indeed observed ||232 |, but 100% to bottomonium 



production. Three basic pQCD production mechanisms, illustrated in figure |IT 
are considered ||233|| : 



Heavy quark fragmentation ||230|| , Z° — > (T)bb. Here the production of a 



primary bb pair is followed by the radiation of a gluon which splits into 
a second bb; a b and b from these two pairs then bind in a colour singlet 
system. It is noteworthy that in addition to the quarkonium two other 
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heavy hadrons occur in this process. 

Gluon fragmentation [231|, Z° — >qqg*, g* — > (bb)gg. The need to remove 



the colour from the bb pair, whilst forming a positive C parity T (or 
Jjijj in the case of cc), requires the emission of two perturbative gluons. 
Thus this is an order-af process and typically the quarkonium state has a 
relatively large transverse momentum. 
• Gluon radiation |234f , Z° — > (bb)gg. Here two gluons are emitted from 
a primary heavy (anti) quark allowing it to recoil and form a colour 
singlet with the other (anti) quark. This results in a very hard, isolated, 
quarkonium state. 

All of the above prompt production mechanisms give comparatively isolated 
quarkonia originating from the interaction point, in contrast to charmonium 
from b decays. Here the bb system is produced in a colour singlet state 
by the emission of perturbative gluons. The dominant process is quark 
fragmentation | 233 |. However the theoretical predictions are significantly low 
compared to the Z° data. Also no evidence for displaced vertices, associated 
with additional heavy quarks, is observed in quarkonium events. A similar 
situation has occurred at the TEVATRON where colour singlet fragmentation 
mechanisms fail to account for the number of observed high-pr quarkonium 
states f23|, |23§ . 

Fortunately recent theoretical developments suggest that colour octet 
quarkonium production may play an important role ||237|| . Here the QQ 
forms a colour octet system from which the colour is leached away (by 
exchange of soft gluons) in an, as yet unspecified, non-perturbative mechanism. 
Two new contributions, also shown in figure [El], arise: gluon fragmentation 
Z° — >qqg*, g* — > (bb); and gluon radiation Z° — »bbg. These new diagrams are 
lower order in a s but suppressed by larger powers of v, the relative velocity of 
the b and b. To calculate these processes requires knowledge of the octet matrix 
elements; until recently these have been taken from fits to the CDF data [[236|| 
though now a potential model based calculation is available [ [23811 . (Ab initio 
lattice calculations have so far only been performed for octet decay matrix 
elements ||239|| and not the technically demanding production matrix elements). 
In the the octet case the dominant fragmentation contribution to T production 
becomes gluon fragmentation |[233j1 . Taking into account all contributions, 
agreement with the Z° data is possible. Although the data are insufficient to 
be able to isolate components due to the individual processes, the lack of an 
observed hard T spectrum does rule out a large contribution from octet gluon 



radiation [|232|| . It is also noteworthy that the octet mechanism predicts a large 



transverse polarization for the vector quarkonium states [p38| , [240|| particularly 
so at high momentum. 
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Colour singlet Colour octet 




heavy quark fragmentation 




gluon fragmentation 




gluon radiation 



Figure 11. The leading production mechanisms for heavy quarkonia 
assuming the colour singlet and octet mechanisms. 
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Looking at the J/ip data in table 4, one can conclude that since b hadron 
decays dominate charmonium production, HERWIG is deficient in this aspect 
of its b decay description. Possible remedies include adding more explicit decay 
channels to the decay tables or refining the partonic b decay model, particularly 
in its treatment of light clusters. Next, because none of the models considered 
include the dominant octet production mechanism the large underestimates 
of T [ |206| , |228| | production are not too surprising. Only ARIADNE, version 
4.09 onwards, been extended to include an approximation to the perturbative 
colour octet production mechanism |241| . This program is therefore the only 
one containing the necessary physics to attempt to describe the T data. 



5-4- Semi-inclusive momentum spectra 

Having discussed the total production rates of identified particles, we now turn 
to their momentum distributions: cr _1 do"/dx. All particle spectra vanish as 
x — > 1 and are expected to vanish as x — > 0, though measurements are rarely 
available at sufficiently low momentum to see indications of this behaviour. 
However changing to the commonly used variable £ = — lnx, a clear 'hump- 
backed plateau' shape (see section |2~5| ) is seen. 



5.4-1- Light hadrons A compilation of many of the measurements of the x 
and £ spectra of identified hadrons at the Z° can be found in |86| and ||181|| . 
References to the more recent measurements are listed in table 6. In addition, 
a recent measurement of the charged particle momentum spectrum is available 
in M 



Table 6. A summary of momentum spectrum measurements that have 





appeared since [86 


]. The letters ADL indicate the contributing experiments. 


Mesons 


Particle 


Reference 


Particle Reference 



V 

u(782) 
r/(958) 
0(1020) 



A 
A 
L 

AL 
D 



|83f 

III 



[189 



[y, igg i 



[190 



£(1520) 
K° 

K*(892)± 
K*(892)° 
K*(1430)° 



D 

A 
A 
D 
D 



&9' 
119' 



Baryons 



Particle 
A 

S(1385) ± 



A 
A 



Reference 

H 

[Hi 



Particle 



S(1530) c 



A 
A 



Reference 

[43] 
[43] 



In section 



2.5.2 



the calculation of the charged particle momentum spectrum 
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was discussed; these results can also be applied to identified hadron spectra |8| . 
Especially in the case of identified hadrons a dependence on the primary quark 
flavour is anticipated: Monte Carlo studies and recent SLD results [j220|| suggest 
that the main effects are at high momentum but that residual effects could also 
occur at low momentum. Consequently investigations have again concentrated 
on the peak position, £*, which occurs at low momentum and for which the 
complete prediction is of the form 

C /• Tin f'% , 

* I V A )\ 2 VA, 

Observe that the cut-off, Qq, dependence only occurs through the first term 
which has the property that F(0) = whilst the scale, Q, dependence only 
arises in the second term: the constant A is the usual (effective) QCD scale. 
In order to create a particle of a given mass it is reasonable to expect that the 
final-state partons will require virtualities of the same magnitude, so that Qq 
is a simple function of the hadron mass ||244|| . Invoking the LPHD concept 
one can hope that these parton level calculations are then sufficient to predict 
an identified hadron's energy spectrum. 

Application of Q24]) has caused some confusion. First, there are small 
differences between experiments in the practical definition of £*. Second, it 
is common to set Qq = A so that the first term vanishes, the limiting spectrum 
case, but then somewhat inconsistently to introduce a new variable Q' Q = A in 
the second term. That is fl2"4"D is replaced by: 




£*=oM£ ( 25 ) 



Given the availability of low energy data it is possible to investigate the In Q 
dependence of £* for the n + , 7r°, K + , K° , p and A hadrons separately [|245| , |243| . 
The data for each hadron species appear to lie on straight lines each of slope 
~ 1/2 but with differing offsets, indicating that F in ( p4f) decreases with 
increasing hadron mass. ALEPH also finds that the linearity of these fits 
can be improved by using JETSET to remove the effects of secondary hadron 
decays ||245|| . For example kaons produced in the weak decays of hard B mesons 
stiffen the kaon spectrum causing a decrease in their £*. However this correction 
procedure is somewhat at odds with the original idea of LPHD which was 
thought to account for such decays. At the Z° it is also possible to study 
the hadron mass dependence of with Q = Mi fixed ||182j , |193| , |243 |, which 
probes the relationship between Qq and m^. Adopting (|25]) OPAL ||182|| claim 
two linear relationships between and (an exponential dependence of Q' 
on one describing the mesons, except pions, and the other the baryons. 
This pattern is well reproduced by JETSET allowing a correction to remove 
secondary decays to be applied after which all the points, including pions, now 
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fall on a single curve. However DELPHI [ 193 | report that mesons and baryons 
lie on two curves corresponding to a negative logarithmic dependence of on 
mi (linear relationship between Q' and wih) which again reduce to a single 
curve after correcting back to primary hadrons using JETSET. 

Whilst the full shape of the momentum spectrum can be calculated within 
pQCD || [63|, [78], [79|] comparison is more commonly made to Monte Carlo 
predictions, particularly so in the case of identified particle momentum spectra, 
reflecting the expected interplay between perturbative and non-perturbative 
contributions. From a practical point of view the primary concern is with a 
correct description of the pion momentum spectra |86|, |182| , |183| , |185| , |245|| which 
dominate the total event multiplicity. These fall by five orders of magnitude 
over the measured momentum range. Here, and in general, the predictions of 
ARIADNE and JETSET are rather similar and typically agree with the data 
a little better than HERWIG, though all suffer problems at very low and very 
high x. Meson spectra are described reasonably well on the whole, though the 
K ± spectrum is notably too soft |]86], |182| , |193| , [245|| which might be attributable 
to an inadequate description of heavy hadron decays. However the description 
for baryons is less satisfactory. In particular the number of fast, predominantly 

In the case 



primary, protons is consistently overestimated ||182j , |193| , 
of ARIADNE and JETSET this may be remedied [43], g|] by invoking a 
suppression of leading uds-baryons , which may have an underlying physics 
motivation 11160 . 



An intriguing possibility is that of an inequality in the fragmentation 
function of strange quarks into protons and antiprotons, D^(x,Q 2 ) ^ 
Dp(x,Q 2 ) [[246 1 . This prediction follows from a possible asymmetry, of non- 
perturbative origin, in the sea quark structure functions of protons ||247|1 and 
application of the 'reciprocity rule' |[248|| . A test of this idea may be possible by 



SLD using a strangeness tag [ 249 1 and exploiting their polarization asymmetry 
to distinguish the quark from the antiquark jet. 



5.4-2. Heavy quark hadrons Hadrons containing heavy (c,b) quarks are 
special since the heavy quark is expected to be principally (if not exclusively) 
of perturbative origin. In practice non-perturbative physics also plays a 
role and the delineation of the two contributions is not clear cut. If the 
fragmentation function can be written as a convolution of perturbative (PT) 
and non-perturbative {NP) parts, D(x) = dp? ® d^p, then it immediately 
follows that (x) = (x)pp ■ (x)np- In leading order pQCD predicts p50|| 
(x) PT = (a s (£ jet )/a s (m Q )) 8 /( 9 ^°); the NLO calculation of d PT (x) J^TJ is 
sharply forward-peaked, vanishing at x — 1. The perturbative result alone is 
too hard {(xb) pt\e=m z /2 ~ 0.8) and a non-perturbative component is required, 
especially in those rare cases (of order 1%) where no gluon radiation occurs at 
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all. 

A simple argument [214] shows that the non-perturbative hadronization 
should also be very hard for heavy quarks due to their inertia, {x)np ~ 
1 — R~ 1 /mQ, where R is a typical hadron size, with the remaining energy 
uniformly distributed in rapidity between ±ln(y / s/2mq). A full expression for 
<Inp{x) depends on the details of the hadronization mechanism assumed, and 
several are available. 



The most commonly adopted standard is due to Peterson et al ||252|| . 



In an independent fragmentation approach, the amplitude for the transition 
(QqO+Q.' i n perturbation theory is proportional to the inverse of the energy 
transfer, assuming a constant matrix element. The fragmentation function is 
then given by the amplitude squared and the appropriate flux factor as: 

D^ X ) = ^(i- 1 --^\ 2 ^Xn^^ iV^VfQ (26) 

X \ X 1 — X J 6q 7T 

Theoretically one predicts eq = R~ 2 /tuq; experimental fits yield small e values 
with e c /e b ~ 10, consistent with (m b /m c ) 2 J25|, |25|, |55], [25|, |57|, |58|. On 



average, (x) = 1 — y/€Q ~ 1 — R^/mq as anticipated above. 

A perceived problem with the Peterson form is its ~ (1— x) 2 behaviour (it is 
sometimes argued that including gluon radiation is equivalent to an additional 
(1 — x) factor) as x — > 1 which is in conflict with the 'reciprocity' rule 



This posits that the x — > 1 behaviour of the Q— >h fragmentation function 
should equal the x — > 1 behaviour of the h— >Q structure function F®(x,Q 2 ). 
Furthermore it is expected from dimensional counting arguments [|259|| that 



as x — > 1, F^(x,Q 2 ) ~ (1 — x) 2ria ~ l where n s is the number of spectators 
in the hadron H. In a meson n s — 1 and therefore the x —* 1 behaviour of 
the fragmentation function should be (1 — x). Two alternatives are available. 
By adopting an explicit meson wavefunction ||260|| and thereby introducing a 
non-trivial matrix element, a 'refined' Peterson form may be derived: 

D%(x)=N(l+x 2 )\ + 4 " 1 i — A^— J -(27) 

\ x 1 — x I V x 1 — xj e 



An earlier approach |261 |, based directly on the structure function analogy, 
leads to: 

D$(x) = (a + I) (a + 2)x a (l - x) (a + l)(a + 2)(1 - x)(28) 

Here a ~ mq so that (x) = 1 — 2/ (a + 3) in accord with expectation. 

In string models two approaches have been elaborated according to whether 
the momentum or space-time aspects are emphasized. In the presence of 
massless quarks, which move along linear light-cone trajectories, these pictures 
are simply related, Ap = nAt, AE = Ap z ; however massive quarks move 
along displaced hyperbolae. In the Lund momentum-space approach the 
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relevant string area is given, as for light quarks, by m\/(xn 2 ) and the 
same reasoning then leads to the LSFF (14). However this gives (x) ~ 
1 — (1 + a)/bm 2 L which is too hard compared to both theoretical expectation 
and experiment ||253| , |254| , [255| , |256| , [257| , |258|| . Adopting a space-time based 
approach | 139 |, a la Artru Mennessier, the string fragments randomly into 
clusters and the relevant string area becomes: 



rn 



2 r 

Q 



2n 2 



m. 



str 



- - 1 - In 









X ) 



(29) 



where m str is the mass of the string segment containing the heavy quark. In 
the Q/Mq — > oo limit this results in a generalization of the LSFF but with no 
a-term present. Elaborating this scheme ||262|| to multiple string breaks along 
a finite size string leads to the following effective fragmentation function: 



1 



X 



l-xV 3 



X 



exp 



-bm 2 ! 

x 



(1-x) 



(30) 



Now m str ~ mi should be identified with the mass of the lightest Q hadron, 
where in fact the cluster mass spectrum peaks. This has a softer spectrum, 
(x) ~ 1 — i? _1 / m Qi an d is very similar to the Peterson form in practice. 

It should be noted that the non-perturbative fragmentation functions 
can be expected to describe effects due to the perturbative emission of soft 
gluons, as found in a parton shower. However effects due to the emission 
of hard gluons cannot be accounted for. In particular when defining x as 
a light-cone momentum fraction this may lead to a difference between the 
reconstructed value x rec , measured along a jet axis, and the primary value x pri , 
generated with respect to the 'string' axis: the Lund model indicates [ 263 1 that 
(x rec -x pri ) b ~ 0.08. 

The experimental measurement of the b fragmentation function is difficult; 
to date it has been based on reconstructed B->DW^(X) decays or the rapidity 
method | 264 |. Also, the interpretation of the results is delicate, and ambiguous 
conclusions have been drawn. Two sources of confusion are the delineation of 
the perturbative and non-perturbative contributions and whether a distinction 
is made between primary and secondary b hadrons cascading down from excited 
states. At the Z° there is substantial production of B* ||200| , |221| , [222|| and the 
four B** [201, 221 , 265] mesons: the primary rate is approximately 1:3:2 for 
B:B*:B**. OPAL ||257|| finds reasonable agreement with the Peterson et al 

and Lund-Bowler 



(p6|), Collins and Spiller (0), Kartvelishvili et al (|2~8| ) and Lund-Bowler ( pCf ) 
fragmentation functions. On the other hand, the ALEPH data [254] would rule 
out the Collins and Spiller form whilst favouring the Kartvelishvili et al form; 
the data also disfavour an untuned version of HERWIG, which is too soft. 
Preliminary studies from DELPHI ||255|| indicate that JETSET with parton 
showers (PS) and the Peterson form, (|26|) , gives the best description of the 
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data but with a peak that is appreciably too wide. Using JETSET PS with 
Lund-Bowler functions, (0), gives a very poor fit as does using JETSET with 
matrix elements and the Peterson function. This last combination gives too 
narrow a peak, indicating the need for a contribution from gluon radiation. 
HERWIG gives a reasonable description of the data but must be tuned to 



avoid too many heavy b clusters. DELPHI ||255|| also give measurements of the 
primary B* and B** fragmenation functions. L3 |256| find agreement within 
errors with the Peterson et al fragmentation function. 



6. Spin phenomena 

As discussed in |L~T| , the primary quarks produced via e + e~— > Z°/7* — > qq 
are highly polarized. Whether, and in what circumstances, this polarization 
survives the hadronization or whether spin-spin forces wash out any memory of 
the initial polarization is an open question. The primary quarks may become 
constituents of unstable baryons or of vector or tensor mesons, the angular 
distribution of whose decay products may be used to extract information about 
spin states. 

In the particular case of leading (large x) spin-1/2 A- type baryons, 
considerable polarization is expected [|266|| , though also see ||267|| . In the 



constituent quark model, the spin of such baryons is carried by the heavy quark, 
with the light diquark system in a spin zero, isospin zero state (though see the 
discussion in section |3.2.5| ). Thus fast Ab particles could carry a substantial 
fraction of the polarization of the initial b quark, with the light diquark system 
carried along as a spectator. Similarly, high-x A c , and A baryons formed in 
fragmentation of s quarks, are expected to be polarized. Since A particles can 
also arise from hadronization of initial u and d quarks, the polarization in this 
case is considerably reduced. 



6.1. Ab polarization 



ALEPH have measured the polarization of Ab baryons ||268|| using semileptonic 
decays, Ab — > /~T7;+charmed hadrons. The method [|269|1 is based on 



measurement of the ratio of the average lepton to average neutrino energy. 



The measured longitudinal polarization, Pa 



-0.23 



+0.25 
-0.21) 



is well below the 



theoretical expectation of —0.69 ± 0.06 although the error is large. This is 
a surprising result and, if confirmed, would indicate the likely existence of 
depolarizing mechanisms in the b quark hadronization. 
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6.2. A polarization 

Because of the parity- violating nature of the decay, A — > pn~ , the distribution 
of the polar angle, 8*, of the proton direction in the A rest frame, relative 
to the A direction in the laboratory, is proportional to 1 + aPi cos 8* , where 
Pl is the degree of longitudinal polarization. The value of the weak decay 
parameter a is well measured . ALEPH ||270|| have exploited this distribution 
to measure a value Pl = —0.32 ± 0.07 for leading (x p > 0.3) A baryons. When 
all sources of A baryons are taken into account, a result of —0.39 ± 0.08 is 
expected if A's containing a primary s quark carry all of its initial polarization. 
Thus the ALEPH measurement is in agreement with standard electroweak 
theory together with the assumption that the initial strange quark polarization 
survives hadronization to become a leading A baryon. 

This result is clearly at odds with the conclusion from the Ab study 
which indicated significant depolarization in the hadronization. Indeed, the 
polarization of the heavier b quark may be expected to survive more easily 
than that of the s quark. 



6.3. Vector meson spin alignment 

Study of spin alignment of vector mesons, particularly at large x, where the 
meson may be expected to contain one of the primary quarks, may provide 
information on the nature of the quark to meson transition. Such analyses 
are normally done in terms of the vector meson helicity density matrix, p\y, 
some of whose elements can be determined by measuring the distribution of 
the vector meson decay products [|271|| ■ The element poo is the fraction of 



mesons which are in the helicity zero state. In strong vector meson decays it is 
not possible to infer separately the values of elements pn and p_i_i, since the 
decay angular distributions are the same for helicity +1 and helicity —1 vector 
mesons. 

In statistical models | [27211 , the fragmentation is assumed to produce extra 
quarks with both helicities equally likely. Parallel alignment of primary and 
secondary quark spins will produce a vector meson with helicity A = ±1. 
If the spins are initially antiparallel, the value of poo will depend on the 
relative probability to produce a vector or a pseudoscalar meson; in this case 
Poo= (1 — P/V)/2, with a maximum value of 1/2 when P/V = and a value 
of 1/3 if there is no suppression of vector mesons. In the model of [ 273 ] 



vector mesons are produced via vector currents q— >qV which conserve the 



quark helicity. The vector meson then has helicity zero. Another model ||274|| . 
for production of leading mesons, assumes multiple emission of soft gluons 
by the fragmenting quark, a process which conserves the quark helicity. The 
leading vector meson may be formed when the leading quark combines with 
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a soft antiquark in a process which results in meson helicity ±1. The basic 
string and cluster models have little to say about vector meson spin alignment, 
although does point out that no alignment may be expected in the simplest 
string picture. 

ALEPH DELPHI f2UD| and OPAL |27§ have measured the alignment 
of B* mesons using the angular distribution of 7 rays from the decay B* — tB'j. 
The measurements agree, and produce a weighted average of poo= 0.33 ± 0.04 
(although they all express the result in terms of a relative contribution of 
longitudinal polarization states). These results, which imply no spin alignment, 
are consistent with simple spin counting and with heavy quark effective theory 
(HQET). They are also in accord with the measurements of the ratio B* to 
B meson production [[200] , 221 , [222|| which imply no suppression of the vector 
state (see |5.2.3| ). 

Results from OPAL [[276|| on lighter vector mesons show deviations from 
Poo= 1/3- For D* 1 * 1 mesons, a value p o = 0.40 ± 0.02 has been measured, 
consistent with lower energy results | 277 |. And for primary 0(1020) mesons at 
x > 0.7, OPAL report an even larger value, poo= 0.54 ± 0.08. These mesons 
therefore appear to be preferentially in the helicity zero state. Measurements of 
some off-diagonal elements of the helicity density matrix show small deviations 
from zero for both mesons. Such non-zero off-diagonal elements are a natural 

m 



consequence of coherence in the hadronization, and are a firm prediction 
of any general model other than independent fragmentation. There is clearly 
more to be learned about hadronization from such measurements. 



7. Correlation phenomena 

Correlations in hadronic systems may be defined as departures from phase 
space in distributions for groups of two or more particles. Such correlations 
may be associated with 

• quantum mechanics — Bose-Einstein and Fermi-Dirac effects 

• hadron dynamics — resonances, reflections and final-state interactions 

• local baryon number conservation - baryon-antibaryon phase space 
correlations 

• local strangeness conservation — strange particle rapidity correlations 

• soft gluon coherence in QCD showers two-particle momentum 
correlations and possibly intermittency 

• fragmentation dynamics - transverse-momentum limitation of phase 
space 

In so far as one wishes to understand the underlying dynamics of the 
hadronization phase, it is important to take into account all of these effects. 
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1.1. Bose-Einstein correlations 

Bose-Einstein correlations (BECs) have been extensively studied in hadronic 
systems from Z° decay. A review may be found in [ f279| ]. The correlations are 
interesting in their own right as a quantum mechanical phenomenon whose 
experimental details can give information on the space-time structure of the 
source of hadrons. The correlations arise from the necessity to symmetrize the 
wavefunction for systems of two or more identical bosons. In most experimental 
analyses, a simple model is assumed where the source of particles is spherical 
with a Gaussian density. Then the two-particle phase space is enhanced by a 
factor C{Q) = 1 + A exp(— Q 2 R 2 ) relative to its density in the absence of the 
correlations. Here Q, the square of the 4- momentum difference between two 
bosons, is the measure of the separation in phase space, A measures the degree 
of coherence in the particle emission (A = corresponds to full coherence) and 
R is the radius of the Gaussian source. 



Table 7. Measurements of Bose-Einstein correlations in Z° decay. 



Particle System R (fm) 



A 



References 



Comments 



K±K± 
KgK§ 



0.65 ±0.16 
0.49 ±0.05 
0.93 ±0.15 
0.62 ±0.05 
0.48 ±0.08 
0.71 ±0.07 
0.55 ±0.14 
0.76 ±0.15 



0.51 ±0.12 
1.06 ±0.17 
0.87 ±0.14 
0.28 ±0.09 
0.82 ±0.27 
1.4 ±0.3 
0.61 ±0.23 
1.14±0.39 



a mj 



D |281 



O [282 



D 1 283 



D 1 284 



A 1 194] 



D 1 284 



O [195 



Direct n 



not corrected for fo 



BECs have been studied at LEP for systems of tt^tc^, tt^h^tc^, tt+tt^tt^, 
K ± K ± and KgKg. A summary of measured values of A and R is given in table 7. 
It is not straightforward to compare the various results, nor to interpret the 
findings. For example, in order to measure the enhancement C(Q) due to BECs 
it is necessary to know the phase space density in their absence, and there are 
several ways to tackle this. In addition particles which arise from resonance or 
weak decays may be removed (if the detector has the capability). The BEC 
effect will be diluted in data samples containing mixtures of different particle 
types, and the purity of pion or kaon samples may be increased, depending 
again on the capabilities of the detector. In the LEP analyses, various different 
approaches to these problems have been used. 

Nevertheless, in all cases it is clear that the model based on a spherical 
source with Gaussian density gives a reasonable fit to the observations, 
although there is no real evidence that other models would not fit equally 
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well. The coherence parameter A varies from close to zero to greater than 
one; the latter is unphysical and possibly indicates deviations from the model. 
The source has a typical size of 1 fm, which at first sight seems inconsistent 
with a picture of hadrons arising from a rapidly expanding linear colour string. 
However the length scale measured by the BECs is not the longitudinal size of 
the string, but the distance in production points for which particles are close 
together in momentum space |15U|| . Recently BEC studies at LEP have been 
extended to three-particle systems | |283| |, and the multiplicity dependence of 
the correlations has been investigated ||285 l. 

Even if all of the measurements had been made in a consistent way, the 
interpretation of the results would not be straightforward: the correlations 
depend in a poorly understood way on final-state interactions, resonance 
production and rescattering of resonance decay products. According to [|153|| 
the situation may be "impossibly complicated" . And ||279|| says "it seems very 
difficult to make progress in studying the Bose-Einstein effect in the context of 
e + e~ physics, and it is not clear to what extent it can be considered a useful 
and interesting activity." 

It is nevertheless important to understand at least the phenomenology 
of the Bose Einstein correlations since they impact on studies of other 
features of hadronic systems. For example, although the effect primarily 
influences systems of identical bosons, in the relatively high-multiplicity, jet- 
like environment of Z° decay, "residual" correlations ||286| arise between pairs 
of unlike particles. The BECs produce a general collimation of the jets and a 
tendency to reduce the mean transverse momentum; this brings all particles 
closer together in momentum space. The effect is to produce a distortion 
of 7r + 7r~ mass spectra, especially at low momentum where the multiplicity 
is highest. This means for example that the use of opposite charge particle 
pairs to determine the phase space in the absence of BECs can result in 
biased values of A and R. An important practical result of the residual 
correlations is the considerable difficulty in inclusive measurements of 7r + 7r~ 
38], |287|| such as the p(770)° and the f (982) whose line shapes 



resonances 



are distorted by the correlations and possibly also by other mechanisms ||286|[ . 

In the process e + e~— > W + W~ — > hadrons at LEP 2, correlations may arise 
between the hadrons from one W and those from the other, an effect recently 



investigated by DELPHI |[288|| . Such an effect could result in a shift of the 
reconstructed W mass in multihadronic W decays ||156|| . A good understanding 
of the role of BECs in hadronization at the Z° may help to reduce uncertainties 
introduced by this effect. 
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7.2. Fermi-Dirac anticorrelations 

While BECs between identical bosons are firmly established experimentally, 
anticorrelations between identical fermions are more difficult to observe. In 
Z° decay, baryons are produced at much lower rates than mesons, and local 
conservation of baryon number suppresses production of identical baryons close 
together in phase space. The only evidence for anticorrelations comes from a 
study by OPAL [[289|| of AA, A A and AA pair production. Close to threshold it 
is found the the spin states of the di-hyperon systems agree with expectations 
of a simple statistical mixture, with no indication of any resonance in the AA 
threshold enhancement. In the case of the identical baryon pairs, there is a 
tentative indication of suppression close to threshold, as would be expected by 
the Pauli exclusion principle. However, further measurements are needed to 
confirm this result. 



7.3. Baryon-antibaryon phase space correlations 



Since baryons, and strange hadrons, are both heavier and less frequently of 
secondary origin than ordinary hadrons investigating their production and 
pairwise correlations appears to offer a more direct probe of the momentum and 
quantum number flow during hadronization. However, allowance must always 
be made for secondary hadrons coming from decays, such as S° — > A — ► p (here 
Hxl on A^p) correlations would be welcome). Baryon number 



further data 162 



conservation implies that a baryon is always accompanied by an antibaryon; 
flavour conservation, via diquark pairs, suggests that the baryons may be 
preferentially particle-antiparticle pairs; and LPHD argues that these baryons 
occur close by in phase space. So far, studies have concentrated on measuring 
the proximity in phase space of p(p) |H |91] , A< A> p|, |09|, g9| and A<p> p2] 



baryon pairs. By introducing an event axis, typically the thrust or sphericity 
axis, correlations can be studied in rapidity, azimuthal angle or polar angle 
with respect to the axis. 

Several models of baryon production are available |[293|| . Section [| discusses 
the independent fragmentation |L6|, |115|| , the diquark [|158|| and popcorn [|133| , 
|160|| (triplet) string options and the cluster ^(J models. In addition a 
possible contibution from direct, 7* — > (qq') (qq')^ diquark production has 
been proposed ||157| , |294| . Also there is the recombination model of hadron 
production |295| . In this approach a spectrum of partons occuring after a 
shower at the fixed scale Ql is convoluted with an explicit wavefunction for a 



particular hadron. In the case of a baryon [|296|1 this gives (c.f. (7)): 
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Db(^Q 2 )= E Jdx 1 dx 2 dx 3 D a aiM (x 1 ,x 2 ,x 3 ,Q 2 ;Ql) 



a 1: a 2 ,a 3 

"ai ,02,03 1 

Here a typical wavefunction is R oc (xiX2X3) 2 5(x — Xi — X 2 — X 3 ), consistent 



X Ra 1 ,a 2 ,a 3 (Xl,X2,X 3 ,X;Ql) (31) 



with reciprocity [ |248|| and quark counting rules | |259| |, and on the assumption 
of uncorrelated partons D a aia2a3 (x u x 2} x 3} Q 2 ]Ql) = Hi D^(xi,Q 2 ;Q 2 ). It 
may be mentioned in passing that the recombination model makes interesting 
predictions for baryon polarizations |297| . 



7.3.1. Rapidity correlations At the Z° strong, short-range correlations in 
rapidity are seen between baryon-antibaryon BB pairs. (Unlike-baryon 
distributions have the like-baryon (BB) distributions subtracted to remove 
secondary correlations due to more than one BB pair in an event.) For 
example, the distribution of baryon rapidity (y) in such events, dn/dy(yB\y^), 
is measured to be compact and centred on y^ JJ3], |162j , |194| , |292| |. In a BB 
pair, given the rapidity of the p (A) there is ~ 70 (50)% probability that 
the p (A) will be found within \y^ — Ub\ < 1(0.6) and vice versa |43|, [292|1 ; a 
marginally weaker short-range Ap correlation is found [|162|| . A much weaker 
long-range anticorrelation is also seen for far forward/backward BB pairs, as 
anticipated from a leading particle effect. In contrast, when BB pairs occur 
in an event, the dn/dy(yB 1 \yB 2 ) distribution is nearly flat with just a weak 
short-range anticorrelation, particularly so away from the central region [[292| 



where phase space constraints become important. Similar, though statistically 
limited, results have been seen at lower energy [|298| , |290 . 



Qualitatively these features are reproduced by both the HERWIG and 
JETSET models, though quantitatively the strength of the correlations is 



overestimated by HERWIG and by JETSET without popcorn |43|, [292|| . The 
recombination model predicts rapidity correlations which are much too weak 
whilst direct diquark production predicts long-range correlations that are too 
strong; both are disfavoured by the data ||209 . 



7.3.2. Azimuthal angle correlations Somewhat weaker correlations are seen 
between the azimuthal angles in baryon pairs at the Z° K3L 194 . BB pairs 



show a tendency towards A0 = 0, though in those pairs which lie out of the 
event plane there is a tendency towards A0 = n. BB pairs show a weaker 
tendency towards A<f) = ir. Two competing mechanisms might be envisaged 
to explain the azimuthal angle correlations. First is a local compensation of 
transverse momentum which leads to an enhancement for A0 = ir. Second is a 
tendency for any off-axis boost to be shared by neighbouring baryons, leading 
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to an enhancement at A(f) = 0. The latter effect might be expected to become 
more important at higher energies where three-jet effects come into play (these 
by definition involve higher transverse momentum scales than hadronization), 
more so for baryons lying in the event plane. At low CoM energy (10 GeV) 



the first, back-to-back, effect appears to dominate [299] but this is seen to 



weaken at higher energy (30 GeV) |[298| , |300| , |301|| and at the Z° the second, 
side- by-side, effect is more important. 

7.3.3. Polar angle correlations Historically a very powerful way to 
discriminate between hadronization models is the orientation of the BB pair 
with respect to the event axis ||301 |. If 8* is the angle between the axis of the BB 



pair and the event axis as seen in the BB pair rest frame then cos 8* is measured 
to be highly forward-backward peaked |43|, |209|| . By utilizing the lepton beam 



polarization information (see section |1 . 1| ) SLD are capable of determining the 
primary quark direction and have recently been able to demonstrate that the 
baryon preferentially follows the quark direction ||291 . 



In string models the colour field typically aligns along the event axis, 
so that provided that the transverse momenta acquired by the baryons are 
small compared to the longitudinal momenta transferred from the string, the 
baryons will retain a strong memory of the string/event axis direction. In 
contrast, because clusters are deemed to be structureless, they have no means 
of retaining any information about an original event axis (they are unpolarized) 
and so decay isotropically in their rest frame. Not surprisingly then the highly 
peaked cos 8* strongly disfavours the essentially flat prediction from HERWIG 
whilst JETSET offers a good description of the data. In HERWIG version 5.7 
an option was introduced to allow non-isotropic decays of clusters containing 
primary quarks; this was implemented to stiffen the momentum spectrum of c 
and b hadrons and has little influence on BB pair polar angle distributions. 

7.4- Strangeness correlations 



In the previous section, [7J|, phase space correlations were discussed for 
baryons; we now turn to similar measurements made on strange hadrons. The 
conservation of strangeness in strong interactions implies that strange hadrons 
are pair-produced during hadronization. The probability that hadron hi is 
accompanied by h 2 is defined as P(hi, h%) = 2 X (n.h lt h 2 )/ ( n h! + n h 2 ) where the 
2 is included because of double counting. Measurements give P(A, A) = 49±6% 
whilst P(A,A) = 13±1% Pirapg^ l, ^(A,K S °) = 17±2% and P(K°, K°) = 



29 ±4% P(S-A) ^40±7%pqp^ and P(~-~") = 4 ± 6% 



Two sources of strange hadrons can be anticipated: leading hadrons 
associated with initial ss quarks and those pair-produced locally in the event's 
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colour field. To help distinguish these possibilities requires information on 
the phase-space correlations between pairs of strange hadrons. A short- 
range rapidity correlation occurs between K°Kg pairs and a slightly weaker 
one between K°A pairs ||194|| indicating a local mechanism for strangeness 
compensation. Also visible are: weaker long-range correlations, as expected 
from leading quarks; and evidence of phase-space suppression, particularly 
when both hadrons are leading. Weak correlations are also seen at A0 = 0, tt 
for centrally produced K°K° and K° A pairs. 

The possible influence of introducing the popcorn baryon production 
mechanism, see section [3.2.5| , is of particular interest for string models. The 
presence of an intermediate meson, BMB tends to soften all correlations. At 
present, taking account of systematic errors, measurements of baryon and 
strangeness correlations are insufficient to place any significant constraint on 
the level of popcorn production required. The BB pair rapidity difference is 
mainly sensitive to the amount of popcorn production, with data favouring 



a substantial component ||209| , [292j| . The substantial rate of S^A pair 
production also favours a high level of popcorn production in order to supply an 
intermediate kaon |190|, [292||. Interestingly the number of B^B^ pairs decreases 



linearly with the amount of popcorn introduced so that their measured 
multiplicity can be (rather simplistically) used to constrain the amount of 
popcorn to around 50% | |194| |. 

A more direct test of the popcorn mechanism is to look at rapidity ordered 
BMB triples ||162|| . Using nearby baryon pairs, |?/b — Ub\ < 1, a probability 
of 7% (25%) is obtained for finding an intervening kaon (pion) in a p'-p) A(A) 
or Ap pair. In addition to there being no enhancement of popcorn-favoured 
AK + p triples, unfavoured AK~p triples are found equally likely. However, in 
events containing a kaon and a close Ap pair, the kaon is found very close in 
rapidity to the baryon pair. 



7. 5. Intermittency 



Intermittency ||302|| , the non-random clustering of particles in phase space, 



is a somewhat obscure phenomenon of uncertain dynamical origin. In 
essence, intermittency corresponds with large, non-statistical, fluctuations in 
the numbers of particles in particular events which are found in narrow rapidity 
bins. It is normally studied by measuring factorial moments of multiplicity 
distributions in rapidity bins. Intermittent behaviour has been observed in 
hadronic systems from Z° decay at both LEP ||303|| and SLC ||304|| as well as in 
lower energy e + e~ collisions. 

The intermittency observed in hadronic Z° decays is in fact reproduced by 
the Lund parton shower Monte Carlo model with string fragmentation. It is 
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a possibility that the self-similarity inherent in the QCD parton shower, with 
its successive q— >qg and g— >gg branchings, is the source of the intermittency 
seen in the distribution of the final-state hadrons. Certainly the JETSET 
model contains no feature explicitly introduced to simulate the dynamics of 
intermittency. 

The appearance of intermittency in the Z° decay data would seem to be 
an ideal opportunity to gain a good understanding of the mechanisms which 
lie behind it. Previously its interpretation in hadropro duct ion experiments has 
been obscured by the complicated nature of the final states and the effects of 
beam and target fragments. 



8. Quark-gluon jet differences 



The determining property of quark and gluon jets is the colour charge of the 
initiator partons, Cp = (N% — l)/(2iV c ) and Ca = N c respectively. Due 
to the gluon's larger charge it should radiate more subsequent gluons in a 
parton shower. This leads to the anticipation that gluon jets, as compared 
to quark jets, will have: a higher multiplicity, softer momentum spectrum 
and wider angular distribution ||305| , |306|| . These global features are largely 
borne out by experimental results ||307|, |72|, [308], |309|, EJTO, 311]. However 



whilst clear differences are now established between gluon and light-quark 
jets, b-jets appear to be rather like gluon jets at Z° energies ||308| , |309| | ; 
the ratios of measured properties typically fall short of the naive asymptotic 
predictions. For example the multiplicity ratio is, in leading order, predicted 
to be C A /Cp = 9/4 TO ; at NNLO this becomes « 2 || |3T| and after 
imposing energy (but not momentum) conservation on the shower this drops 
to ~ 1.6 — 1.8 ||313|| (the exact predictions depend on jet energy and the scale 
used for a s ). The measured ratio, which is seen to be sensitive to the precise 
jet definition, is typically in the range 1.1 — 1.3, though OPAL has obtained 
1.55 in an event hemisphere-based analysis |310|| . It should be borne in mind 
that the these basic predictions are made for back-to-back pairs of quark or 
gluon jets whilst growing evidence suggests that the relative topology of a jet 
is important in determining the appropriate scale ||3 1 1|| . 

The above discussion relates to the perturbative properties of quark 
and gluon jets. Current cluster and string models of hadronization make 
no distinction between whether a set of final-state partons arose from a 
fragmenting primary quark or gluon; they are treated the same. This does 
allow the possibility of a 'leading particle' effect (for example one should 
expect more leading kaons in an s- quark jet than a gluon jet) but no other 

. An interesting possibility is that some isoscalar 
may contain a significant gg component and hence 



'anomalous' effects 
states, r], r]', 
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might appear more frequently as leading particles in gluon jets. These particles 
are very often primary hadrons coming directly from a cluster or string. At 
present the possibility of hadrons containing gluons is not allowed for and only 
quark constituents are considered. Indeed in the cluster framework gluons are 
split into light qq pairs whilst in the string approach they represent energy- 
momentum 'kinks' on a string. However modifying the models to accommodate 
gluonium would not, a priori, appear to pose significant problems of principle 
or practice. 

A more radical scenario is offered by the independent fragmentation model 
of Peterson and Walsh [314] . This is based on the suggestion that a gluon 
is attached to an octet colour flux tube and quarks to triplet flux tubes (see 
section |3.2|) . A gluon now fragments into a sequence of isoscalar gg or gqqg 
clusters leading to a prediction of greatly enhanced rj, rj' etc production and 
harder momentum spectra as compared to quark jets. 

A third alternative scenario for gluon jet hadronization is provided within 



the recombination model ||295|| discussed in section [773]. Calculations predict: 
a softer pion spectrum |pl5|| , an enhanced multiplicity ratio (n r? /)/(n 7r o) [ [316| 



and enhanced baryon production ||296|| in gluon jets compared to quark jets. 

The L3 Collaboration have reported tentative indications that rj production 
is enhanced in gluon jets [|186|| . Studying the lowest energy — gluon — jet in 
three-jet events they see a a harder 77 momentum spectrum than predicted by 
both the HERWIG and JETSET Monte Carlos with an enhancement in the 
ratio (n v ) / (n n o). The use of a ratio takes into account the established increase 
in multiplicity found in gluon compared to quark jets and is designed to make 
the measurement sensitive to any additional enhancement or suppression. The 
Monte Carlos provide a satisfactory description of the spectrum and {n v ) / (n w o) 
ratio in quark jets. DELPHI have studied the production rates of kaons, A 
and S in multi-jet events [ 2X2 1 . They find that the relative yields of strange 
hadrons in multi-jet, normalized to two-jet, events is constant in events with 
widely separated jets but favours increased production in multi-jet events at 
small resolutions, particularly so for kaons. In a more direct study DELPHI 
have looked at identified particle production rates in actual quark (natural 
flavour mix) and gluon jets, normalized to the charged multiplicity in the jet: 
the double ratios ((%) / (n c h)) g / '((^h) / '{n c h)) g I [31711 - The evidence suggests 
values for the double ratios of approximately 1.1 for the K°, 0.9 for K + , 1.2 
for p and 1.4 for A. In general these results are in qualitative agreement with 
the Monte Carlos but quantitatively the deviations from unity are larger than 
typically predicted, apart for the K° result where a small suppression was 
expected. However the errors are relatively large. OPAL have also reported 
preliminary studies of the double ratios ||318 |. For K° and mesons they report 
a slight, < 10%, increase in the relative production rates in gluon jets, whilst for 
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p and A baryons they measure a significant 30-40% increase. These enhanced 
production rates are presumably at the expense of pion production. At present 
more measurements are required before hard conclusions can be drawn. 

The possibility of enhanced identified hadron production in gluon as 
compared to quark jets has received previous attention. The results for 
mesons, in particular hard r\ relative to ir°, production are not conclusive. At 
y/s = 10 GeV ARGUS [319] saw no evidence for enhanced r\ or production in 
the continuum, 7* — >qq. Also comparison has been made between continuum 
and T(l, 25') — >ggg events; here Crystal Ball [|320|| see no enhancement whilst 
DASP-II, CLEO and ARGUS ||321|| see a slight enhancement for a number 
of mesons. At \fs = 30 GeV JADE ||322|| reported very weak (statistically 
insignificant) evidence for a small enhancement in the (n v )/(n n o) ratio in 
acollinear, gluon rich, events. The situation is clearer for baryons. The 
production rate of baryons in T(l, 25*) decay (gluon dominated) to continuum 
(quark dominated) events shows an excess of 200-300% [P21|| . 

In the context of conventional Monte Carlos a number of partial 
explanations have been offered. A study using JETSET indicates that 
the relative production rates of mesons in quark and gluon jets is energy 
independent and just less than one, whilst for baryons it is 20-25% larger 
and shows a slight increase with jet energy |[318|| . This latter effect may 
be attributed to an edge effect associated with the suppression of leading 
baryons. In cluster models it has been argued that in gluon rich environments, 
the topology, rather than any intrinsic properties of the jets, leads to heavier 
clusters and hence larger baryon production rates \'5'2'i\. At the T it has also 
been emphasized [[324 1 that secondary decays are important, and that 40% of 
the events are in the continuum to b and c quarks (see figure [l]) which, being 
heavy, 'eat up' the available phase space for baryon production. 



9. Outlook 



That Monte Carlo event generators, solidly based on sound physics, are 
essential in modern high energy experiments, from detector conception to data 
analysis, ought not to be forgotten. It therefore almost goes without saying 
that the reliability of Monte Carlo predictions should be a prime concern if 
only for mere practical reasons. To this end the large event rates and pristine 
conditions available in hadronic Z° decays at e + e~ colliders play a particularly 
important role. Here, as nowhere else, the physics assumption built into the 
Monte Carlo models can be confronted with ever more exacting tests. The 
relatively complex conditions associated with initial-state hadrons have largely 
precluded this activity using ep and pp data, thereby making physicists reliant 
on the quality of the Z° data. Two caveats to the wider application of the 
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models are the questions of the reliability of the factorization theorem in 
perturbative QCD and the presumption of universality in the hadronization 
process. 

Focusing on the models for the hadronization processes one may ask what 
is left that should be done with the Z° data and what might be learned from it. 
We list a number of topics, perhaps not all of which can be done with available 
data (but one can always hope for more Z° data): 

• A study of the transverse momentum distributions of identified particles, 



especially pions where low mass effects and correlations may occur [|147 
- In string models the fact that the predicted width of the Gaussian 



p_i_-distribution, y/t/vr, proves too narrow has raised questions about the 
adequacy of a tunnelling mechanism explanation: is unresolved (non- 
perturbative) gluon emission [|145|| the real explanation? 
An attempt to establish properties of the relatively rare, directly produced 
pions — Pions are special particles by virtue of their nature as Goldstone 
bosons (of the chiral symmetry); this mandates them to have small masses 
which in a string model implies a very small size, in fact uniquely less than 
a string's width. 

Further measurements of orbitally excited mesons and baryons — Are the 
production rates for these states too high to be compatible with the string 
model? 

A search for D-wave states, such as the Kg (1780) - These could 
help elucidate the role of the wavefunction in determining a hadron's 
production rate. 

A study of f (980) production as a function of rapidity and of total event 
multiplicity — This would test the Gribov confinement scenario ||215| . 
A search for deuterium production — This has been reported previously 
in studies at the T [|325|| and is expected, on the basis of a string model 
calculation ||326|| , at the level of 5 x 10~ 5 per hadronic Z° event. 
Further measurements of strangeness suppression (7 S ) which ideally can 
be directly compared to those available in ep collisions, particularly those 



associated with the current region of the Breit frame |88[ — If it turns out 



that jg e 7^ 7| p then effort should go to establishing any other differences. 
A measurement of the s-quark to p and p fragmentation functions — Is 
there a measurable difference as suggested ||246|| by the possible asymmetry 
in sea quark structure functions? 

Measurements of identified particle production rates in quark and gluon 
jets — Are there measurable differences and, if so, what is the mechanism? 
A search for glueball candidates, particularly in the gluon-rich 
environment provided by gluon jets. 
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Further measurements of leading baryon polarization - Is there a 
depolarizing mechanism for the b baryons which does not apply to those 
arising from s quark fragmentation? 

Measurements of the helicity density matrix for the light and heavy vector 
mesons — Does hadronization produce spin-aligned vector mesons, and 
does any alignment depend on x7 Are some off-diagonal elements non- 
zero as predicted by coherent hadronization? Spin physics has a long 
history of producing surprises. 

More detailed studies of intermittency — Can it be firmly established that 
the phenomenon is due to the self-similar evolution of the initial qq state 
via a QCD parton shower? 

A search for direct evidence of popcorn-type baryon production via, for 
example, further study of AK + p type correlations — Three-body cluster 
decays would also induce such correlations but perhaps with different 
intensity 

A search for Q~p type correlations — These are possible in a generalized 
popcorn mechanism although they will be hard to find experimentally 
because of low rates. 

An attempt to establish directly the existence of the 'dead cone' [ |110|| 
in non-leading particles, by removing the leading particles using fully 
reconstructed b-hadron decays. 

Direct measurement of the rates and momentum spectra of hadrons 
produced in b-quark events to ensure adequate descriptions of the weak 
decays in the Monte Carlos — b hadrons contribute a tenth of all particles 
in hadronic Z° decays, and significantly more at large x. 
A tuning of the Monte Carlo models to Z° data, incorporating colour 
rearrangement - - This is important in order to determine the effects 
of colour rearrangement and to provide constraints which may prove 
significant for later W mass measurements. 

A simultaneous tuning of the Monte Carlo programs to the Z° data and 
lower energy data, particularly from PEP and PETRA — At lower energies 
the relative contribution to event properties made by hadronization is more 
important, whilst the influence of the perturbative shower can be tested 
using the CoM energy dependence of observables. 

An 'ultimate' tuning of the Monte Carlo models to the final Z° data — This 
will be an invaluable service to experimentalists and theorists allowing the 
experience gained from LEP1/SLC to be applied at future machines. 
A continuation of the search for tests which discriminate between 
the competing models of the non-perturbative physics underlying the 
hadronization, and continuing development of these models. 
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In summary, the e + e~— > Z°/7* — > hadrons data have already provided 
a wealth of information on the phenomena of parton hadronization. Of the 
available models of the non-perturbative physics involved, the Lund string 
model, as implemented in JETSET, has met with most success, particularly in 
the baryon sector, and most notably in its prediction of the angular distribution 
of correlated baryon-antibaryon pairs in their rest frame. However, the model 
has many free parameters and consequently has little predictive power. But 
the parameters are not arbitrary — most are based on incomplete knowledge 
of physics. Therefore it could be argued that their values, when fully tuned to 
reproduce observations, provide important information about hadronic physics. 
On the other hand, while the weight of evidence tends to favour the string 
picture, the other models, particularly the cluster model of HERWIG, are not 
dead, and more analyses of the existing and future data are essential to provide 
further discrimination between the models and to help elucidate the physics of 
hadronization. 
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